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Abstract — Interference alignment has emerged as a powerful 
tool in the analysis of multi-user networks. Despite considerable 
recent progress, the capacity region of the Gaussian A'-user 
interference channel is still unknown in general, in part due 
to the challenges associated with alignment on the signal scale 
using lattice codes. This paper develops a new framework for 
lattice interference alignment, based on the compute-and-forward 
approach. Within this framework, each receiver decodes by first 
recovering two or more linear combinations of the transmitted 
codewords with integer-valued coefficients and then solving these 
equations for its desired codeword. For the special case of 
symmetric channel gains, this framework is used to derive the 
approximate sum capacity of the Gaussian interference channel, 
up to an explicitly defined outage set of the channel gains. The 
key contributions are the capacity lower bounds for the weak 
through strong interference regimes, where each receiver should 
jointly decode its own codeword along with part of the interfering 
codewords. As part of the analysis, it is shown that decoding 
K linear combinations of the codewords can approach the sum 
capacity of the it-user Gaussian multiple-access channel up to a 
gap of no more than ^ log K bits. 



I. Introduction 

Handling interference efficiently is a major challenge in 
multi-user wireless communication. Recently, it has become 
clear that this challenge can sometimes be overcome via 
interference alignment fT), Q. For instance, consider the 
if -user Gaussian interference channel, where K transmitter- 
receiver pairs wish to communicate simultaneously. Through 
the use of clever encoding strategies, it is possible to align 
the transmitted signals so that each receiver only observes 
its desired signal along with a single effective interferer. As 
a result, each user can achieve roughly half the rate that 
would be available were there no interference whatsoever, 
i.e., K/2 degrees-of-freedom (DoF) are available. However, 
many schemes, such as the Cadambe-Jafar framework [2] and 
ergodic interference alignment A3], require a large number 
of independent channel realizations to achieve near-perfect 
alignment. In certain settings, this level of channel diversity 
may not be attainable; ideally, we would like to achieve 
alignment over a single channel realization. 

The capacity region of the (static) Gaussian i\"-user interfer- 
ence channel [4] is unknown in general, although significant 
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progress has been made recently, in part due to the discovery 
of interference alignment and the shift from exact capacity 
results to capacity approximations ||5|-||7). It has been shown 
by Motahari et al. that K/2 degrees-of-freedom are achievable 
for almost all channel realizations [8 1 but it is an open question 
as to whether this result translates to real gains outside of the 
very high signal-to-noise ratio (SNR) regime. One promising 
direction is the use of lattice codes ll9l- lfTTI . as they can 
enable alignment on the signal scale. By taking advantage of 
the fact that the sum of lattice codewords is itself a lattice 
codeword, a receiver can treat several users as one effective 
user, thereby reducing the number of effective interferers. A 
compelling example of this approach is the derivation of the 
approximate capacity of the many-to-one interference channel 
by Bresler, Parekh, and Tse [7|. For fully connected channels, 
much less is known, owing to the difficulty of choosing lattices 
that simultaneously align at several receivers. 

In some cases, focusing on the special case of symmetric 
channel gains has yielded important insights. For instance, 
in the two-user case, Etkin, Tse, and Wang [5| used the 
symmetric interference channel to develop the notion of gener- 
alized degrees-of-freedom. This in turn revealed five operating 
regimes, based on relative interference strength: 

• Noisy: Each receiver treats interference as noise, which 
is optimal for sufficiently weak interference 11121 - 11141 . 

• Weak and Moderately Weak: Each transmitter sends a 
public and a private codeword following the scheme of 
Han and Kobayashi fl5l . Each receiver jointly decodes 
both public codewords and its desired private codeword 
while treating the interfering private codeword as noise. 

• Strong: Each receiver jointly decodes both users' code- 
words. This regime and its capacity was discovered by 
Sato IT6ll as well as Han and Kobayashi fl5l . 

• Very Strong: Each receiver decodes and subtracts the 
interference before recovering its desired codeword. This 
regime and its capacity was discovered by Carleial [17]. 

Using these regimes as a guideline, they were able to approx- 
imate the capacity region to within half a bit. 

In this paper, we focus on the special case of the symmetric 
Gaussian K-user interference channel. Each receiver observes 



yfc = x fc 



9 x " 



Zfe 



where is the codeword sent by the kth transmitter, g is 
the cross-channel gain, and z^. is additive white Gaussian 
noise. If each transmitter uses the same lattice codebook, 
then each receiver observes an effective two-user multiple- 
access channel (MAC). The first effective transmitter sends the 
desired codeword and the second sends the sum of interfering 
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codewords, which is itself a codeword. A natural approach is 
to first decode and cancel this sum of codewords J2 m ^k Xm ' 
leaving only the desired signal plus noise. In the very strong 
interference regime, this approach is optimal lfj"8l but it does 
not suffice in general. In our setting, joint decoding must cope 
with the fact that all users employ the same lattice codebook, 
and distinct pairs of lattice codewords may result in the same 
sum. 

Our main contribution is the derivation of new achievable 
rate regions for the weak, moderately weak, and strong inter- 
ference regimes, for which no capacity approaching schemes 
were known in the literature. This is enabled using a new tech- 
nique, based on compute-and-forward lfl9l . that allows each 
receiver to jointly decode its desired codeword and the sum of 
the interfering codewords. We also propose a generalization of 
the two-user Han-Kobayashi scheme fT3Tl in which each user 
transmits a superposition of one private lattice codeword and 
one public lattice codeword. Each receiver decodes its private 
and public lattice codewords as well as linear combinations 
of all interfering public lattice codewords while treating the 
interfering private lattice codewords as noise. Using these new 
tools, we can mimic the achievable scheme of [5], and obtain 
an approximation of the sum capacity in all regimes. In the 
weak interference regime, the approximation is valid for all 
channel gains, whereas in the strong and moderately weak 
interference regimes, it is valid for all channel gains except 
for an outage set. The measure of this outage set can be made 
as small as desired, but this comes at the expense of increasing 
the gap between the inner and outer bounds on the capacity. 
Our capacity approximation closely follows the basic shape of 
the generalized degrees-of-freedom (shown in Figure [3], which 
was derived for this channel by Jafar and Vishwanath [20 1 . 
Their approach is inspired by the deterministic model [6 1, and 
relies upon coding across an infinite number of signal levels. 
As a result, their coding scheme cannot be directly translated 
into a finite SNR result, as existing techniques for decoding 
the sum of codewords incur a penalty for each additional 
codeword layer [19|. 

In its original incarnation, compute-and-forward makes it 
possible for relays in a network to decode integer combina- 
tions of the transmitted codewords and send them towards a 
destination fI91 - The effective SNR at each relay is determined 
by how closely the integer coefficients match the channel 
coefficients. If together the relays recover a full-rank set of 
linear combinations, the destination can recover all of the 
messages. In our setting, each receiver decodes two linear 
combinations of the transmitted codewords, 



anx fc 
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where an, 012,021, and 022 are integer-valued coefficients. If 
the vectors ai = [an ai2\ T and &2 = [021 a22] T are linearly 
independent, then each receiver can solve for its desired code- 
word Xfc. Through this approach, we can derive closed-form 
lower bounds on the performance of joint decoding. A direct 
analysis of joint decoding, although possible [21], presents 
additional difficulties that in turn lead to looser bounds. 
As part of the derivation of our sum capacity bounds, we 



develop a new decoding framework, the compute-and-forward 
transform, that may be of independent interest. Consider a 
Gaussian J^-user MAC where the encoders employ lattice 
codes. Under our framework, the receiver attempts to recover 
K linearly independent equations of the transmitted codewords 
and then solve them for its desired messages. This transforms 
the MAC into an effective multiple-input multiple-output 
(MIMO) channel with an integer-valued channel matrix. A 
striking phenomenon we discover is that while the computation 
rate for each of these K equations is very sensitive to the exact 
channel gains, the sum of their computation rates is equal to 
the sum capacity of the MAC up to a constant gap, independent 
of the channel gains and the SNR (see Figure [TJ. To give 
this an operational meaning, the computation rate for each 
equation is associated with a single transmitter's codeword. 
The receiver successively decodes the equations, ordered by 
decreasing computation rate. As we will show, each decoded 
equation can be used to remove its associated codeword from 
the receiver's observation in a way that does not alter the 
effective SNR for the remaining equations. 
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Fig. 1. Computation rates for the best two linearly independent equations 
vs. h for the channel y = xi + /1X2 + z at SNR=40dB. The sum of these 
computation rates is nearly equal to the multiple-access sum capacity. All 
rates are normalized by this sum capacity l/21og(l + (1 + h 2 )SNR). 



A. Related Work 

Interference alignment has generated a great deal of excite- 
ment, due to the promise of higher throughputs in wireless 
networks |fl~), 10 as well as other applications, including 
coding for distributed storage ||221 . See the recent monograph 
by Jafar for a comprehensive survey ll23l . Of particular note is 
a series of recent papers that delineate the degrees-of-freedom 
limits of linear beamforming strategies for alignment over a 
finite number of channel realizations [24|, [25 1. Beamforming 
strategies can only approach perfect alignment asymptotically, 
whereas lattice-based schemes can achieve K/2 degrees-of- 
freedom over a single channel realization [8 1. However, lattice- 
based alignment at finite SNR has to date been limited to 
special cases, such as symmetric |18|, l2ll . Il26j . integer [27|, 
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and many-to-one interference channels [7], [28]. Capacity ap- 
proximations are also available for one-to-many [7 1 and cyclic 
interference channels ||29l , although these coding schemes do 
not employ alignment. 

The results in this paper are connected to the recent work 
of Niesen and Maddah-Ali Il301 . which proposes a finer- 
grain deterministic model for interference networks that is 
rich enough to capture the phenomenon of alignment. Via this 
model, they derive the approximate capacity region of the two- 
user Gaussian X-channel. As in our work, they characterize 
the performance of joint decoding over a MAC where all 
users employ the same linear code. While both approaches 
yield "constant-gap" approximations, our framework naturally 
yields achievable rate expressions that are easy to plot and 
often much closer to the upper bound than the constant-gap 
analysis suggests. 

Bandemer and El Gamal have recently proposed a class of 
three-user deterministic channels where the interfering signals 
are passed through a function on their way to the receiver, 
which, in a certain sense, models interference alignment [ 3 1 1 . 
They develop a new rate region based on interference decoding 
for this model. In a recent paper, Wu, Shamai, and Verdu have 
derived a general formula for the degrees-of-freedom of the 
K-user Gaussian interference channel via Renyi's information 
dimension l32ll . 

Nested lattice codes have been thoroughly studied as a 
framework for efficient source and channel coding with side 
information J9), iflOl . [33]. Recently, it has become clear 
that the inherent linear structure of lattices can enable many 
interesting new schemes, including distributed dirty paper 
coding [34 1, distributed source coding of linear functions [35|- 
(37), distributed antenna systems [38 1-[40|, and physical-layer 
network coding JT9), BTl - ll44l . to name a few. The origins 
of these schemes can be traced to the work of Korner and 
Marton B31 . who showed that linear binning is optimal for the 
distributed compression of the parity of a doubly symmetric 
binary source. 

The remainder of the paper is organized as follows. Section 
|n]gives a formal problem statement as well as the approximate 
sum capacity of the Gaussian K-user interference channel. 
Next, SectionHIIlprovides a brief review of nested lattice codes 
and the compute-and-forward strategy. Section [IV] shows how 
the Gaussian multiple-access sum capacity can be approached 
within a constant gap using compute-and-forward combined 
with a type of successive interference cancellation. It also 
shows how an effective multiple-access channel emerges in 
the context of lattice interference alignment. The upper bounds 
needed to establish our approximate sum capacity result are 
reviewed in Section [V] Afterwards, Section |Vl] develops our 
two achievable schemes and develops closed-form expressions 
for each interference regime. Finally, Section fVTTI derives the 
degrees-of-freedom attained by our strategy. 

II. Symmetric Gaussian JT-User Interference 
Channel 

A. Problem Statement 

We begin with some notational conventions. We will denote 
column vectors with boldface lowercase letters and matrices 



with boldface uppercase letters. Fo r instance, a e T, K and 

A e Z Kx K . Let ||a|| = \Jj2k=i a l denote the £ 2 -norm of the 
vector a. Also, let denote the zero vector and Ikxk denote 
the identity matrix of size K. We use [•] to denote rounding to 
the nearest integer, |_-J to denote the floor operation and [•] for 
the ceiling operation. In general, the letters a and b are used 
in this paper whenever the variables they describe are integer 
valued. All logarithms are to base 2. We also occasionally use 
the notation log + (x) = max(0, log(a;)). 
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Fig. 2. Block diagram of a symmetric Gaussian K-user interference channel. 

There are K transmitter-receiver pairs that wish to simul- 
taneously communicate across a shared channel over n time 
slots, where the channel gains are constant over all n channel 
uses. We assume a real-valued channel model throughout. 

Definition 1 (Messages): Each transmitter has a 
message Wk drawn independently and uniformly over 
{1,2,..., 2™^™}. 

Definition 2 (Encoders): Each transmitter is equipped with 
an encoder, £ k : {1, 2, . . . , 2 nRsYM } -> 1", that maps its 
message into a length-n channel input = Sf.(wk) that 
satisfies the power constraint, 

||xfc|| 2 < nSNR 

where SNR > is the signal-to-noise ratio. 

Definition 3 (Channel Model): The channel output at each 
receiver is a noisy linear combination of its desired signal and 
the sum of the interfering terms, of the form 



Yk = x fc 



9 x„ 



Zfc , 



(1) 



where g > parametrizes the interference strength and Zfc 
is an i.i.d. Gaussian vector with mean and variance 1. We 
define the interference-to-noise ratio to be 

INR = ff 2 SNR 
and the interference level to be 

a log(INR) 



a 



log(SNR) 



Remark 1: Note that our definition of INR ignores the fact 
that there are K — 1 interferers observed at each receiver. 
This is for two reasons. First, this definition parallels that of 
the two-user case [5], which will make it easier to compare 
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the two rate regions. Second, the receivers will often be able 
to treat the interference as stemming from a single effective 
transmitter, via interference alignment. 

Definition 4 (Decoders): Each receiver is equipped with a 
decoder, V k : R" ->■ {1, 2, . . . , 2 nRs ™}, that produces an 
estimate Wk = T>k(yk) of its desired message Wk- 

Definition 5 (Symmetric Capacity): A symmetric rate 
i?SYM is achievable if, for any e > and n large enough, 
there exist encoders and decoders that can attain probability 
of error at most e, 



Pr ({wi ^ wi} U • • • U {w K ^ W K }\ 



< e 



The symmetric capacity Csym is the supremum of all achiev- 
able symmetric rates. 

Remark 2: Due to the symmetry of the channel, the sym- 
metric capacity is equal to the sum capacity, normalized by 
the number of users. To see this, assume that the users employ 
different rates and that a rate tuple (i?i, i?2, ■ • ■ , Rk) is 
achievable. Since each transmitter-receiver pair sees the same 
effective channel, we can simply exchange the encoders and 
decoders to achieve the rate tuple (ii^nN, R w (2), ■ ■ ■ ,Rir(K)) 
for any permutation tt. By time-sharing across all permu- 
tations, we find that each user can achieve J2k=i ^ fc > 
corresponding to a symmetric rate. Thus, the sum of any 
achievable rate tuple is upper bounded by if Csym- 

Definition 6 (Generalized Degrees-of-Freedom): The gen- 
eralized degrees-of-freedom (GDoF) specifies the fraction of 
the point-to-point Gaussian capacity that can be attained per 
user for a given interference level a > as SNR tends to 
infinity, 

d(a) = lim , SYM . 

V ' SNR-^oo | log(l + SNR) 

B. Approximate Sum Capacity 

As shown by Jafar and Vishwanath ||20l Theorem 3.1], the 
GDoF of the symmetric if -user interference channel is iden- 
tical to that of the two-user channel, except for a singularity 

at a = 1, 



d(a) 



1- 

a 

1- 
j_ 

K 
a 
2 

1 



< a < | 



§ < a < 1 
a = 1 

1 < a < 2 (strong) 
a > 2 (very strong). 



(noisy) 
(weak) 

(moderately weak) 



See Figure [3] for a plot. Notice that since SNR is taken to 
infinity, the GDoF characterization treats all channel gains g 
that do not scale with SNR as a single point at a = 1. A finer 
view of this regime is possible at high SNR by simply setting 
g to be some fixed value and then taking SNR to infinity, 
corresponding to the standard notion of degrees-of-freedom. 
Surprisingly, this degrees-of-freedom characterization is ev- 
erywhere discontinuous with respect to g [46 1. This presents 
an obstacle towards a clean capacity approximation at finite 
SNR. 




Fig. 3. Generalized degrees-of-freedom for the symmetric Gaussian K -user 
interference channel. 



To overcome this difficulty, our approximations allow for the 
possibility of an outage set, which is explicitly characterized. 
Specifically, in the regime around a = 1, our capacity results 
take the following shape: for any constant c > 0, the capacity 
is approximated within at most c + 9 + log if bits over 
the entire range of SNR, and all channel gains g, except 
for a set of measure ^t(c) which vanishes rapidly with c. 
This type of capacity approximation has also been used by 
Niesen and Maddah-Ali for the two-user Gaussian X channel 
ll30ll and seems to arise from the capacity region itself, not 
just the lower bound. That is, it appears that the capacity 
may in fact simultaneously vary rapidly with the fine scale 
of the channel gains (e.g., the distance to an appropriately 
scaled integer) and slowly on the coarse scale (e.g., relative 
interference strength). In the high SNR limit, this behavior 
shows up as a discontinuity on the rationals but, at reasonable 
SNRs, our achievable scheme shows that this variation is in 
fact fairly smooth. The theorem below captures our capacity 
approximations in a simple form. 

Theorem 1: The symmetric capacity of the symmetric 
Gaussian if -user interference channel can be lower and upper 
bounded as follows: 

• Noisy Interference Regime, < a < ^, 



- loE 

2 S 



1 



1 + INR 



1 



< Csym < - log 



1 



SNR 
1 + INR 



1 



Weak Interference Regime, i < a < |, 

i log + (INR) -I- log(if ) < C SY m < \ log+(INR) + 1 
for all channel gains. 

Moderately Weak Interference Regime, I < a < 1, 



1 



log" 1 



SNR 



c - 8 - log(if ) 
< Csym < ^ log" 1 



SNR x 
v/lNR, 
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(c) SNR = 50dB (d) SNR = 65dB 

Fig. 4. Upper and lower bounds on the sum capacity of a 3-user symmetric Gaussian interference channel with respect to the cross-gain g. The upper bound 
(red line) is given by j421 and the lower bound (black line) is the maximum of the achievable rates from Theorem [9] and Corollary [5] which were computed 
numerically, and Theorem 1101 The lower bounds from Theorem [TJ are not plotted in this figure. For reference, we have also plotted the rate achievable via 
time-division (dotted blue line). 



for all channel gains except for an outage set of measure 
fi < 2~ c for any c > 0. 

• Strong Interference Regime, 1 < a < 2, 

i log+(INR) - | - 3 < C7 SYM < \ log+(INR) + 1 

for all channel gains except for an outage set whose mea- 
sure is a fraction of 2~ c of the interval 1 < \g\ < vSNR, 
for any c > 0. 

• Very Strong Interference Regime, a > 2, 

~ log(l + SNR) - 1 < C SYM < ~ log(l + SNR) 



III. Preliminaries 

In this section we give some basic definitions and results 
that will be extensively used in the sequel. 



A. K-user Gaussian MAC 

Consider the A'-user Gaussian MAC 

K 

y = hkXk + z ' ( 2 ) 

it=i 

where the vector h = [hi ■ ■ ■ hx\ T G i K represents the 
channel gains, x^, 6 W\ k = 1, . . . , K, are the channel inputs, 
z € K n is additive white Gaussian noise (AWGN) with zero 
mean and unit variance and y e M. n is the channel output. 
Without loss of generality, we assume all K users are subject 
to the same power constrain^] 

||x fc || 2 < nSNR, k = l,...,K. (3) 

The capacity region of the channel (O is known (see e.g., Il47l 
Theorem 15.3.6]) to be the set of all rate tuples [R\, . . . , Rk) 

'As otherwise the different powers can be absorbed into the channel gains. 
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satisfying 



^i? fe <ilog[l + SNR^ 

kes \ kes 



(4) 



for all subsets S C {1,...,K}. The achievability part 
of the capacity theorem is established using i.i.d. Gaussian 
codebooks for all users. Motivated by lattice interference 
alignment, we are interested in establishing the achievability 
of certain rate tuples under the constraint that the codebooks 
employed by the K users form a chain of nested lattice codes. 

B. Nested Lattice Codes 

We employ the nested lattice framework originally proposed 
in iflOl . A lattice A is a discrete subgroup of M™ which is 
closed under reflection and real addition. Formally, for any 
ti,t 2 G A, we have that — ti, -t 2 G A and ti +t 2 G A. Note 
that by definition the zero vector is always a member of the 
lattice. Any lattice A in K" is spanned by some n x n matrix 
G such that 

A = {t = Gq : q G Z n }. 

We say that a lattice is full-rank if its spanning matrix G is 
full-rank. 

We denote the nearest neighbor quantizer associated with 
the lattice A by 



Qa(x) = argmin ||x - t| 



(5) 



The Voronoi region of A, denoted by V, is the set of all 
points in M" which are quantized to the zero vector, where 
ties in © are broken in a systematic manner. The modulo 
operation returns the quantization error w.r.t. the lattice, 

[x] mod A = x — <2a(x), 
and satisfies the distributive law, 

[a[x] mod A + b[y] mod A] mod A = [ax + by] mod A, 

for all a, b G Z. 

A lattice A is said to be nested in Ai if A C Ai. The coding 
schemes presented in this paper utilize a chain of K + l nested 
lattices satisfying 



A C A K C • •• C Ai. 



(6) 



From these lattices, we construct K codebooks, one for each 
user. Specifically, user k is allocated the codebook Ck = 
AgriA n V, where V is the Voronoi region of A and the function 
9(k) : {1, ...,K} — > {1, ...,K} maps between users and 
lattices. The rate of each codebook Ck is 

R k = -loglA^A,) nvl. 

n 

User k encodes its message into a lattice point from its 
codebook, t/. G Ck- Each user also has a randorql dither vector 
dfc which is generated independently and uniformly over V. 

2 It can be shown that these random dithers can be replaced with determin- 
istic ones, meaning that no common randomness is required. 



These dithers are made available to the decoder. The signal 
transmitted by user k is 

Xfc = [tfc - dfc] mod A. 

Remark 3: The nested lattice construction from iflCfl em- 
ploys Construction A. To create each fine lattice, this proce- 
dure first embeds codewords drawn from a linear code into 
the unit cube, and then applies the generator matrix for the 
coarse lattice A. As shown in ifTOI . this ensemble of nested 
lattice codes can approach the capacity of a point-to-point 
Gaussian channel. If the integers Z" are selected as the coarse 
lattice, the resulting nested lattice code is equivalent to a 
linear code coupled with a pulse amplitude modulation (PAM) 
constellation. Furthermore, the mod A operation simplifies to 
the quantization error from rounding to the integers. It can be 
shown that the cost of this simplification is only the shaping 
gain, which corresponds to at most 1/2 log(27re/12) ~ 0.255 
bits per channel use. 

C. Compute-and-Forward 





Xl 


Si 






V —*■ v 



mod A 











£k 





Fig. 5. Compute-and-forward on a Gaussian multiple-access channel. 
The transmitters send lattice points t k and the receiver decodes an integer 
combination of them, modulo the coarse lattice A. The rate is determined by 
how closely the equation coefficients a k match the channel coefficients h k . 

Our objective is to communicate over the MAC using 
the compute-and-forward scheme from |fl9l . To this end, the 
receiver first decodes a set of K lattice equations with linearly 
independent coefficient vectors. Afterwards, it solves this set 
of equations for the transmitted lattice points. Assume the 
receiver is interested in decoding the lattice equation 



U-l 



mod A 



with coefficient vector a = [oi ••• ax] G Z . Following 
the scheme of [19], the receiver scales the observation y by a 
factor f3, removes the dithers, and reduces modulo A to get 



K 



Py + ^ akdk 



k=l 



mod A 



K 



K 



K 



afcXfc + ^ °fe d fe + J~]{Phk - a,t)x fe + f3z 



k=l 



k=l 



k=l 



[v + z eff (h, a, f})] mod A, 



mod A 
(7) 
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where 



K 



z eff (h, a, (3) = y^Qgfofc - a k )x k + j3z 



(8) 



k=l 



is effective noise. From 1(191 , we have that z e ff(h, a, j3) is 
statistically independent of v and its effective variance, defined 
as 

a e 2 ff (h,a,/?)^iE||z eff (h,a,/3)|| 2 (9) 

n 

is 

a^(h, a, P) = \\0h- a|| 2 ■ SNR + /3 2 . (10) 

Let k* — minfe :afc ^o be the index of the densest lattice 
participating in the lattice equation v. The receiver produces an 
estimate for v by applying to s the lattice quantizer associated 
with Afc», 



v = [Q\ h » (s)] mod A. 



(11) 



Let Vfc» be the Voronoi region of A&», and define the error 
probability 



Pr(v^v) <Pr(z eff (h,a,/3) ^ V fc .; 



(12) 



The next theorem summarizes and reformulates relevant re- 
sults from Sections IV.C, IV.D, and V.A of fl9l . 

Theorem 2: For any e > and n large enough there exists 
a chain of n-dimensional nested lattices A C A# C • • • C 
Ai forming the set of codebooks C%, . . . ,Ck having rates 
Ri, . . . , Rk and satisfying the power constraint ((3} such that: 
(a) For all channel vectors h £ H K and coefficient vectors 
a £ Z K , the average error probability in decoding the 



lattice equation v 



mod A of transmitted 



lattice points t/. £ C k can be made smaller than e so long 
as the message rates do not exceed the computation rate, 



Rk < i? CO mp(h, a,/3) 



1 / SNR 

2 ° g U 2 ff (h,a,/3) 



(13) 



for all k such that a k ^ and some /3 £ K. 

(b) The codebooks Cx, . . . , Ck are isomorphic to some set 
of linear codebooks Ci, . . . ,Ck over the finite field Z p , 
where p is a sufficiently large prime number. 

(c) For the same p, the equation [p ■ t] mod A = holds 
Vt £ A fe , k= 1,...,K. 

Corollary 1: Given K lattice equations V = [vi • • • Vk] 
with coefficient vectors A = [ai • • • &k] T , the lattice points 



ti, 



. , tx can be recovered if [A] mod p is full-rank over 



Remark 4: Note that it is also possible to map both the 
messages and the lattice equations into an appropriately chosen 
finite field. That is, the messages can be written as vectors with 
elements that take values in a prime-sized finite field, and the 
equations are linear combinations of the messages over the 
same finite field. See |fl9l for more details. 

It follows from Theorem l2l|a) that in order to maximize the 
computation rate i? comp (h, a, (3) for a given coefficient vector, 
one has to minimize c 2 ff (h, a, (3) over /3. It is seen from (TTOb 
that the expression for cr 2 ff (h, a, /3) is equal to the mean 



squared error (MSE) for linear estimation of X = 2~2k=i a kX k 
from Y — Y^,k=i hkX k + Z where {X k }^ =1 are i.i.d. 
random variables with zero mean and variance SNR and Z 
is statistically independent of {Xk)^_ t with zero mean and 
unit variance. Hence the minimizing value of j3 is the linear 
minimum mean squared error (MMSE) estimation coefficient 
of X from Y. This value of j3 was found in |fl9l Theorem 2] 
and the MSE is given by 

^eff( h > a ) - min cr 2 ff (h, a, /3) 



SNR 

SNR a 7 f l /x 



SNRfh^a 



T„\2 



l + SNR||h|| 2 , 
SNR hh T 
l + SNR||h|| 



= a T (SNR- 1 I Kx K + hh T ) 'a 



(SNR-^^xK 



hh a 



(14) 
(15) 



where (fl4l can be verified using Woodbury's matrix identity 
(i.e., the Matrix Inversion Lemma) 11481 Thm 18.2.8]. Accord- 
ingly, we define 



-Rcom P (h, a) = maxfl, 



comp(h, a, p) 

SNR 



2 l0g V^ 2 ff (h,a) 



(16) 



The following definition identifies the K linearly independent 
coefficient vectors which yield the highest computation rates. 

Definition 7: We say that an ordered set of linearly inde- 
pendent integer coefficient vectors {a 1; . . . ,&k} with corre- 
sponding computation rates i? CO mp,fe = -Rcomp(h, a^-) is optimal 



if R, 



comp,l 



> • • • > i? C omp,if and for any k — 1, . . . , K and any 



set of integer coefficient vectors {ai, . . . , a*} of rank k 



min i?com P (h,a£) < i? comPi fc, 



or equivalently 

^max fc cr e 2 ff (h,a £ ) > cr e 2 ff (h,a fc ). 

Note that this set is not unique. For example, if {ai, . . . , a^} 
is an optimal set of coefficient vectors, so is the set 
{-ai, . . . ,-a K }- 

Remark 5: Several recent papers have proposed families 
of constellations and codes that are well-suited for low- 
complexity implementations of compute-and-forward [39 1, 
EH, Il49l - ll52"l . These codes could serve as building blocks 
for a practical implementation of our alignment scheme. 



D. Numerical Evaluations 

The optimal coefficient vectors and computation rates from 
Definition [7] play an important role in the achievable rate re- 
gions derived in this paper. The problem of determining the op- 
timal coefficient vectors is that of finding the set of K linearly 
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independent integer-valued vectors that minimizes the effec- 
tive noise (fl~5T >. This problem is equivalent to finding the short- 
est K linearly independent lattice vectors in the lattice^ A(F) 

spanned by the matrix F = (SNR _1 Ik x k + hh T ) 1/2 . 

It is shown in |fl9l Lemma 1] that only integer vectors a £ 
% K that satisfy the condition 



< 1 



h|| 2 SNR 



(17) 



yield positive rates. Therefore, in our considerations it suffices 
to enumerate all integer vectors (other than the zero vector) 
that satisfy ([17l l, and then exhaustively search over these 
vectors in order to find the optimal set. At moderate values of 
SNR this task is computationally reasonable. Nevertheless, it 
is sometimes simpler to find a set of short linearly independent 
lattice vectors in A(F), which is not necessarily optimal, in 
order to obtain lower bounds on the set of optimal computation 
rates. A simple low-complexity algorithm for computing a 
short lattice basis (which forms a set of K linearly independent 
lattice vectors) is the LLL algorithm [53 J In producing the 
figures for this paper we have employed the LLL algorithm, 
meaning that the plotted achievable rates in Figure |4] are in 
fact lower bounds on the rates given by Theorems [9] and Q~T] 
We note that a similar procedure for finding the optimal 
coefficient vectors was also described in P4l . where the 
optimal coefficient vectors are termed dominated solutions. 

IV. Multiple- Access via Compute- and-Forward 

This section introduces a new coding technique for reliable 
communication over the if-user Gaussian multiple-access 
channel. This technique is based on the receiver decoding K 
linearly independent equations of the transmitted codewords, 
and then solving them for obtaining the messages transmitted 
by each user. We begin this section with a high-level overview 
of the scheme, which is illustrated in Figures [6] and [7] 

Each user k maps its message to a lattice point tfc in its 
codebook Ck and transmits a dithered version of it. The K 
lattice codebooks utilized by the different users form a chain 
of nested lattices as in ©. Assume for now that the users 
are ordered with descending rates R\ > R2 > • • • > Rk, i.e., 
6(k) — fc for k — 1, . . . , K. The receiver, which sees a noisy 
real-valued linear combination of the transmitted codewords, 
begins by decoding the integer-valued linear combination 
vi = ai m t m ] mod A which yields the highest computa- 
tion rate i? comp .i. Using the compute-and-forward framework, 
this is possible if Ri < i? comp ,i. Then, it proceeds to decode 
the equation v 2 = a 2m t m ] mod A which yields the second 
highest computation rate i? C omp,2- In general, ti participates in 
this equation and the condition for correct decoding of v 2 is 
therefore i?i < i? C omp,2- Nevertheless, this condition can be 
relaxed using the first equation vi that was already decoded. 
Specifically, after appropriate scaling of the channel's output 
and dithers removal, the receiver has a noisy observation 

S2 = [V2 + z e ff(h, a 2 )] mod A 

3 Notice that this A"-dimensional lattice is induced by the channel matrix, 
not the n-dimensional coding scheme. 

4 Pseudocode for the LLL algorithm can be found, e.g., in 1541 . 



of the desired equation v 2 . If ti participates in vi, it 
is possible to cancel out ti from the second equation 
by adding a scaled version of vi to s 2 . Namely, the re- 
ceiver adds r 2 iVi to s 2 , where r 2 i is an integer chosen 
such that [(an + r 2 ia 2 i)] mod p = 0, which assures that 
[(an + r 2 ia 2 i)ti] mod A = for any ti G C\. After 
reducing modA this yields 

sf = [v 2 + r 2 iv a + z eff (h, a 2 )] mod A 
= [v a + z eff (h, a 2 )] mod A, 

where ti does not participate in v 2 . Since the effective noise 
z e ff(h, a 2 ) is unchanged by this process, the receiver can 
decode v 2 as long as i? 2 < ^comp.2- Now, the receiver 
can obtain v 2 by subtracting r 2 ivi from v 2 and reducing 
mod A0 The receiver decodes the remaining equations in a 
similar manner, i.e., before decoding the fcth equation Vfc with 
computation rate i? CO mp.fc the receiver adds to 

Sfc = [v fc + z e ff(h, a fe )] mod A 

an integer-valued linear combination X^=i r ki^i mod A of 
the lattice equations that were already decoded. The coeffi- 
cients in the linear combination are chosen such that the effect 
of ti, . . . , tfe_i is canceled out from Vfc. Assuming that such 
coefficients {r^i, . . . , ^-1} exist, the receiver can decode 
Vfc = Vfc + Y^tZl r k£v e mod A as long as R k < R com p,k- 

Lemma 12 stated in the sequel, establishes that for any set 
of K linearly independent coefficient vectors {ai, . . . , a^} 
there indeed always exist integer- valued coefficients {r^} such 
that in each fcth decoding step the receiver can cancel out 
k — 1 lattice points from the desired equation Vfc, using the 
previously decoded equations {vi, . . . , Vfc_i}. The procedure 
for finding these coefficients is reminiscent of the Gaussian 
elimination procedure of a full-rank matrix. One of the basic 
operations in Gaussian elimination is row switching. In our 
considerations, this would correspond to using an equation 
that was not decoded yet for eliminating lattice points from 
another equation. Since our successive cancelation procedure 
only uses decoded equations, this is not possible. Therefore, 
a major difference between our procedure for finding a good 
set of coefficients {rij} and Gaussian elimination is that row 
switching is not permitted. This incurs a constraint on the order 
in which we cancel out users from equations. Nevertheless, 
there always exists at least one order of successive cancelation 
that is possible. In other words, we can always cancel out the 
effect of k — 1 users from Vfc using the decoded equation 
{vi, . . . , Vfc_i}, but we cannot always control which of the 
K users to cancel. As a result, there always exists at least 
one permutation vector tt such that all K equations can be 
decoded as long as 



R-n{k) < Rcomp.ki k — 1, 



(18) 



It follows that a sum-rate of J\=i i? C omp,fc is achievable over 
the .KT-user MAC with our scheme, where all users are using 

5 The operation of extracting V2 from V2 is in fact not necessary as the 
receiver is only interested in decoding any set of K linearly independent 
equations. We describe this step only to simplify the exposition of the scheme. 
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nested lattice codebooks. As we shall see, this sum rate is 
within a constant gap, smaller than K /2\og(K) bits, from the 
sum capacity of the MAC, for any channel gains and SNR. 

A. The Compute-and-Forward Transform 

We first introduce a transformation of a MAC to a multiple- 
input multiple-output (MIMO) mod-A channel, where the 
K x K channel matrix is integer-valued. This transforma- 
tion, dubbed the compute-and-forward transform, will play an 
important role in our decoding scheme for the interference 
channel. 

Definition 8: Let {ai, . . . , 8lk } be a set of optimal integer 
coefficient vectors (see Definition [7J, /3i , . . . , Pk the cor- 
responding optimal scaling factors, and i? CO mp,i > • • • > 
-Rcomp,A' the corresponding optimal computation rates. We 
define the compute-and-forward transform of the MAC with 
nested lattice codes as 



S = 



si 




where we have written the channel output y, dithers d^, and 
lattice codewords tk as length-n row vectors. We also denote 



[ai 



&k] T and Z eff = [z^ 



Remark 6: The transform is not unique as the set of optimal 
integer coefficient vectors is not unique. Nevertheless, the set 
of optimal computation rates is unique. As we shall see, the set 
of optimal computation rates dictates the rates attained over the 
transformed channel. Therefore, we use the term the compute- 
and-forward transform of the channel, with the understanding 
that although there may be multiple options for the transform, 
they are all equivalent. 

The fcth output S/. of the transformed channel corresponds 
to a lattice equation plus effective noise. Due to Theorem [2] 
each such lattice equation can be reliably decoded as long 
as all lattice points participating in it belong to codes Ck of 
rates smaller than i? com p,fc- We now lower bound the sum of 
K optimal computation rates, and in the sequel we show that 
this sum can be translated to a valid MAC sum rate. 

Theorem 3: The sum of optimal computation rates is lower 
bounded by 

f>c 0mp , fc > \ log (1 + ||h|| 2 SNR) - j \og(K) . (20) 
fe=i 

We will need the following definition for the proof. Note that 
the lattice referred to below is over the K dimensions induced 
by the channel vector h g K , rather than the n dimensions 
used for coding. 



Definition 9 (Successive minima): Let A(F) be a full-rank 
lattice in M A ' spanned by the matrix F G R KxK . For fc = 
1, . . . , K, we define the fcth successive minimum as 

A fe (F) =inf jr : dim (span (a(F) f)B(0, r))) > k\ 

where 13(0, r) = {xet Jf : ||x|| < r} is the closed ball of 
radius r around 0. In words, the fcth successive minimum of 
a lattice is the minimal radius of a ball centered around that 
contains fc linearly independent lattice points. 

The product of successive minima can be upper bounded us- 
ing the following well-known theorem due to Minkowski lf55l 
Theorem 1.5]. 

Theorem 4 (Minkowski): For any lattice A(F) which is 
spanned by a full-rank K x K matrix F 

K 

nA£(F)<A^|det(F)| 2 . (21) 

k=l 

We are now ready to prove Theorem [3] 

Proof of Theorem \3] Let A(F) be a lattice spanned 

by the matrix F = (SNR _1 I KxK + hh T )~ 1/2 , and 
let Ai(F), . . . , Aif(F) be its K successive minima. Let 
ai,...,aif denote the optimal coefficient vectors. By Defi- 
nition|7] Definition |9] and (O we have ||F a fc || = A fc (F) for 
fc = 1, . . . , K. The sum of optimal computation rates is 

K 

^ ^comp(h, afc) 



K 

^ ^ -^comp.fc 
k=l 



k=l 
K 



/c=i 

K 
~2 



SNR 



cr e 2 ff (h,a fc ) 



log (SNR) — i log (jJjF^II 2 ) 
= |log(SNR)-ilog^J]A 2 (F)y 



Applying Theorem [4] to the product JIaLi ^10^) yi e Ws 
K K 

Ei?com P ,fe > y log (SNR) - -log(^|det(F)| 2 ) . (22) 

k=l 

Using Sylvester's determinant identity (see e.g., [48 1) 
det(I^ xK + SNR hh T ) = det(l + ||h|| 2 SNR), 
we have that 

SNR A 



|det(F)r 



(23) 



1 + ||h|| 2 SNR ■ 

Substituting d23l into (l22l proves the theorem. ■ 
Remark 7: It is possible to avoid the loss of the con- 
stant factor K /2 log K in (f20b using successive compute-and- 
forward, as described in [1561 . However, in this case the 
operational interpretation of the sum of computation rates 
becomes more complicated than that described in the sequel. 
We will elaborate on this issue in future work. 

Next, we give an operational meaning to the K optimal 
computation rates. 
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Wi- 



W 2 - 



-di 
ti 1 

t 2 



mod A 



C 



W K - 




mod A 


x 2 hi \ 






h K 


mod A 








E a Jfmdm Y, r KmV m ~J2 r Km^ 

►(+) — 4— I^Mt^ 



mod A 



mod A 


tl 








- 




mod A 


t 2 


/r 1 










mod A 


i K 


K 





-Wl 



■W2 



Fig. 6. System diagram of the nested lattice encoding and decoding operations employed as part of the compute-and-forward transform. Each message 
is mapped to a lattice codeword according to codebook C^, dithered, and transmitted as x^. The multiple-access channel scales codeword k by /ij. and 
outputs the sum plus Gaussian noise z. The decoder attempts to recover K linearly independent equations with coefficients A = {afe m }. For the figure, we 
have assumed that R\ > R2 > ■ • • > Rk and that R^ < i?comp(h, a^, f}^). To decode the first equation vi = [J^ ii m t m ] mod A, the receiver scales y 
by 0i, removes the dithers, quantizes using Qa-l > an d takes mod A. For the second equation V2 = E a 2mt m ] mod A, the decoder scales by /?2, removes 
the dithers, and then eliminates the lattice point ti using the first equation vi so that the rate of the remaining lattice points is at most R^. It then quantizes 
using Qa 2 > a dds back in vi, and takes mod A. Decoding proceeds in this fashion, using a form of successive interference cancellation to keep the rates 
of the lattice points below the computation rates. Afterwards, the receiver solves for the original lattice points by multiplying by A -1 and taking mod A. 
Finally, it maps these estimates of the transmitted lattice points back to the corresponding messages. 



101- 



W2- 



C 2 



w K - 



C 



K 



mod A 



mod A 



mod A 



z eff (h,ai,/3i) 

^4- — > 

z e ff(h,a 2 ,/3 2 ) 



QAi 




mod A 





Vl 




z e ff(h,aK, Pk) J2 - rKmV 

-^(+) k+Mqa^U©— I 



mod A 



A- 



mod A 


tl 


zr 1 










mod A 


t 2 


z - 2 










mod A 


i K 


'-K 





■Wl 



~W2 



Fig. 7. Effective MIMO channel induced by the compute-and-forward transform of a Gaussian multiple-access channel. The channel output y = ^fc x fc + z 
is converted into K linearly independent lattice equations = Ecifc m t m ] mod A plus effective noise z e ff(h, , /3j.) = /3^z + — a fe) x fc- As 

in Figure |6] these lattice equations can be decoded using a version of successive cancellation. 



B. Multiple-Access Sum Capacity to within a Constant Gap 

We now show that the compute-and-forward transform can 
be used for achieving several rate tuples within a constant 
gap from the boundary of the capacity region of the A'-user 
MAC. The main technique used for establishing this result is 
successive cancelation. Namely, each decoded lattice equation 
is used to cancel out the effect of one user from the lattice 
equations that have yet to be decoded. We first illustrate the 
coding scheme by an example, and then formalize our result 
in Theorem 

Example 1: Consider the two-user MAC 

y = V5xi + x 2 + z, 

at SNR = 15dB. It can be shown using < TT3T > and dl6l ) that the 
compute-and-forward transform of this channel is 

si \ \( 2 1 \ / ti \ / z eff ,i 
s 2 ) \ 3 1 ) I t 2 ) + I z eff , 2 



mod A 



with i? comp .i ~ 2.409 bits and i? comp . 2 ~ 1.372 bits. Note that 

(i?comp,i+-Rcom P , 2 )/(l/21og(l + ||h|| 2 SNR)) ~ 0.998. We use 



a chain of three nested lattices A C A 2 C Ai that satisfy the 
conditions of Theorem |2] in order to construct the codebooks 
C\ = Ai n V with rate R\ arbitrarily close to i? C omp,i for user 
1 and £ 2 = A 2 n V with rate R 2 arbitrarily close to i? C omp, 2 
for user 2. 



From Theorem |2jtaJ, we know that vi = [2t% + 1 2 ] mod A 
can be decoded from Si since i?i and i? 2 are smaller than 
^comp.i- However, Theorem |2] does not guarantee that the 
second lattice equation v 2 = [3ti +t 2 ] mod A can be decoded 
from s 2 since the first user employs a codebook with a rate 
Ri ps i? C omp,i which is higher than the second computation 
rate i? C omp, 2 - To circumvent this issue, we use the decoded 
equation Vi as side information in order to cancel out the 
lattice point ti € Ai n V from s 2 . Note that Theorem |2]|cJ 
guarantees that [p • tfc] mod A = 0, k = 1,2 for some 
sufficiently large prime number p. Let 2 _1 e Z be an integer 
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that satisfies [2 1 • 2] modp = 1. The receiver computes 

s 2 = [s 2 - 3 ■ 2 _1 vi] mod A 

= [(3 - 3 • 2" 1 • 2)ti + (1 - 3 • 2~ 1 )t 2 + z e ff, 2 ] mod A 
== [[1 — 3 • 2 _1 ] mod p ■ t 2 + z e ff, 2 ] mod A 



= [ai2 • t 2 + Zeff.a] mod A, 



(24) 



where 612 = [1 — 3-2 -1 ] mod p. Step (a) in d24i i follows from 
the distributive law. Step (b) follows since 3 — 3-2 _1 -2 = M-p 
for some M £ Z. Thus, 

[(3 - 3 • 2 _1 • 2)ti] mod A = [M-p - t x ] mod A 

= [A/ • [p ■ ti] mod A] mod A 
= 0, 

where the last equality is justified by Theorem 121c}. 

Now only t 2 participates in the lattice equation v 2 = 
[fii 2 t 2 ] mod A and, since i? 2 is smaller than -R C omp.2, The- 
orem [2] guarantees that it can be decoded from s|'. This is 
accomplished by quantizing onto A 2 and reducing modulo A, 

v 2 = [QA^sf)] mod A. 

After decoding both lattice equations vi and v 2 the receiver 
can solve for the transmitted lattice points ti and t 2 , as the two 
equations are full -rank over Z p . We have therefore shown that 
the rate region i?i < i? C om P ,i and i? 2 < i? CO mp,2 is achievable. 
In a similar manner, we can show that the rate region R\ < 
^comp.2 and i? 2 < i?comp.i is achievable with this scheme. 

In order to formally characterize the achievable rate region, 
we will need the following definition which identifies the 
orders for which successive cancelation can be performed. 

Definition 10: For a full-rank K x K matrix A with 
integer-valued entries we define the pseudo-triangularization 
process, which transforms the matrix A to a matrix 
A which is upper triangular up to column permutation 
7r = [7r(l) 7r(2) ••■ tt(K)]. This is accomplished by left- 
multiplying A by a lower triangular matrix L with unit 
diagonal, such that A = LA is upper triangular up to column 
permutation tt. Although the matrix A is integer valued, the 
matrices L and A need not necessarily be integer valued. 
Note that the pseudo-triangularization process is reminiscent 
of Gaussian elimination except that row switching and row 
multiplication are prohibited. 

Example 2: The 2x2 matrix 



A = 



2 1 

3 1 



from Example [TJ can be pseudo-triangularized with two differ- 
ent permutation vectors 



1 



-5 



[1 2], 



1 
-1 1 

2 1 
1 



TT = \2 11 



Remark 8: Any full-rank matrix can be triangularized using 
the Gaussian elimination process, and therefore any full-rank 
matrix can be pseudo-triangularized with at least one permu- 
tation vector tt. In particular, since for any MAC the integer- 
valued matrix A from the compute-and-forward transform is 
full-rank, it can always be pseudo-triangularized with at least 
one permutation vector tt. There are full-rank matrices that can 
be pseudo-triangularized with several different permutation 
vectors, such as A from Example [2] However, there are also 
full-rank matrices A that can be pseudo-triangularized with 
only one permutation vector tt. An example of such a matrix 
is the identity matrix Ikxk- 

The next theorem gives an achievable rate region for the 
MAC under the compute-and-forward transform. The proof 
is given in Appendix [A] and follows along the same lines as 
Example [TJ . 

Theorem 5: Consider the MAC ©. For any e > and n 
large enough, there exists a chain of n-dimensional nested 
lattices A C h K C ••• C A; forming the set of codebooks 
£1, . . . , Ck with rates R\, . . . , Rk such that for all h € M. K , 
if: 

1) each user k encodes its message using the codebook Ck, 

2) the integer-valued matrix from the compute-and-forward 
transform of the MAC © can be pseudo-triangularized 
with the permutation vector tt, and the optimal computa- 
tion rates are -R CO m P ,i > • • • > Rcom P ,K, 

3) all rates Ri, . . . , Rk satisfy 



Rk < R 



comp,7r 1 (k) 



, for k = 1 , . . . , K 



(25) 



where 7r _1 is the inverse permutation vector of tt, 
then all messages can be decoded with error probability 
smaller than e. 

Combining Theorems [3] and [5] gives the following theorem. 

Theorem 6: The sum rate achieved by the compute-and- 
forward transform has a gap of no more than K /2 log K bits 
from the sum capacity of the MAC. 

Proof: Let i? C omp,i > • • • > ^comp.A' be the optimal 
computation rates in the compute-and-forward transform of 
the MAC ©. The integer- valued matrix from the compute- 
and-forward transform can be pseudo-triangularized with at 
least one permutation vector tt. By Theorem [5] the rate tuple 



Ri- — R, 



comp,7r 1 (k) 



— S, for k = 1 , . . . , K 



(26) 



is achievable for any <5 > 0. For this rate tuple we have 



A" 



fe=i k=i 

K 

k=l 



comp,7r -1 (fc) 



5) 



comp,fc 



KS 



> I log (1 + ||h|| 2 SNR) 



ylog(*Q 



K6, (27) 
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where (l27T i follows from Theorem [3] Since this is true for any 
S > 0, the result follows. ■ 

C. Effective Multiple-Access Channel 

A channel that often arises in the context of lattice interfer- 
ence alignment is a A-user Gaussian multiple-access channel 
(MAC) with integer-valued ratios between some of the channel 
coefficients. Specifically, the output of such a channel can be 
written as 



1=1 \meK( I 



(28) 



where K\, . . . , JCl are disjoint subsets of {1, . . . , K}. We 
assume that the b m G Z are non-zero integers, which opens 
up the possibility of lattice alignment. 

The channel d28l > may describe the signal seen by a receiver 
in an interference network, perhaps after appropriate precoding 
at the transmitters. In such networks, each receiver is only 
interested in the messages from some of the users while 
the others act as interferers. Hence, it is beneficial to align 
several interfering users into one effective interferer, by taking 
advantage of the fact that the sum of lattice codewords is itself 
a lattice codeword. 

Definition 11 (Effective users): For the MAC specified 
by ( |28l >, we define L effective users 



x e ff . 



= 1, 



• L. 



Definition 12 (Effective MAC): The if -user MAC (Eg) in- 
duces the effective L-user MAC 



L 

1=1 



(29) 



with the vector of effective channel coefficients g = 
[<?i ' ' ' 9l] T <= K L - The effective channel is further char- 
acterized by the effective users' weights 



h 2 



= E ^ 

for £ = 1, . . . , L, and the effective (diagonal) weight matrix 
B^diag(6 e 2 ffil ,...,6 e 2 ffiL ). (30) 

Definition 13 (Effective lattice points): Let t m be the lat- 
tice point transmitted by user m. We define the effective lattice 
point corresponding to effective user I as 



bmXr 



mod A. 



Let 9 e fi(£) = min me jc £ 0(m) (where #(•) is the mapping 
between users and fine lattices defined in Section iHlT i be the 
index of the densest lattice contributing to t e ff/. Since all 
lattices are nested, it follows that t e ff.£ 6 A^^). 

Example 3 (Symmetric K-user interference channel): 
Consider the symmetric A-user interference channel ([T}. 
The channel seen by the fcth receiver is of the form of 



with gi = 1, (? 2 = 9, K-i = {k}, = {!,... ,K} \ k and 
b m = 1 Vm = 1, . . . , K. If each of the K users transmits a 
single codeword drawn from a common nested lattice code, 
the channel becomes an effective two-user MAC as 

Vfc = X eff ,fci + (?X e ff,fc2 + Zfc, 

where the effective users are Xeff.fci = X& and ~x. e s,k2 = 
2~2m^k x m' anc ^ me effective users' weights are b 2 s x = 1 and 
b 2 S2 = K — 1. The effective lattice points are t e fr : ki — tfc and 
t e ff,fe2 = E m ^ tm] mod A. 

Our achievable schemes for the symmetric K-user in- 
terference channel, developed in Section [VI] are based on 
transforming the A-user MAC seen by each receiver into an 
effective MAC with less effective users. We will develop two 
schemes: One transforms the channel into an effective two- 
user MAC as in the example above. The other, which mimics 
the Han-Kobayashi approach, transforms the channel into an 
effective three-user MAC. 

When lattice interference alignment schemes are designed 
properly, the message intended for the receiver is mapped into 
a unique effective user, while multiple interfering users are 
folded into a smaller number of effective users. In this case, it 
suffices for the receiver to decode only the L effective lattice 
points corresponding to the effective users, rather than the K 
lattice points transmitted by all users. In our considerations, 
the effective lattice points are recovered by first decoding L 
lattice equations of the type 



1=1 meKi 
L 



mod A 



teff.< 



.1=1 



mod A 



(31) 



with linearly independent coefficient vectors, and then solving 

for trff,!, . . . , t e ff,L. 

As in Section [Till in order to decode a lattice equation v, 
the receiver first scales its observation by a factor /?, removes 
the dithers and reduces modulo A which yields 



/3y + at 2J b m d m mod A 

1=1 m£K t 
L L 

^ a£X eff! /j + 2J a t 2_j bmd 

L (=1 t=l rnelCt 

L 

+ J~](figi - a£)x eff ,£ + /3z 



mod A 



= [v + z eff (g, a, (i, {b m })} mod A, 

where 

L 

z e ff (g, a, P, {b m }) = f3y - 2J atx eff . e 



(32) 



^2{f3ge - a e ) &m x m + Z (33) 

1=1 melCe 
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is effective noise which is statistically independent of v. Its 
effective variance is 

L 

a 2 eS (g, a, 8, B) = SNR - a e ) 2 b 2 fV + 8 2 , (34) 

l=\ 

where B is defined in (f30b . Let £* = mm£- ae ^Q 8 e ff(£) be 
the index of the densest lattice participating in the lattice 
equation v. Since all lattices are nested, then v G At*. The 
receiver produces an estimate for v by applying to s the lattice 
quantizer associated with Ai», 

v = [Q Ai , (s)] mod A. (35) 

It follows from Theorem [2] that there exists a chain of K + 1 
nested lattices which allows to decode v with a vanishing error 
probability so long as 



R m < i?comp(g, a, 8, B) 



loe 



SNR 



<7^(g,a,/3,B) 



(36) 



for all me \J e:a ^ K, e . 

The expression for o~ 2 ff(g, a, /3, B) is equal to the MSE for 
linear estimation of X — 2~2e=i a t%-i from Y = £« =1 9iXe+ 
Z where {Xe}f =1 are statistically independent random vari- 
ables with zero mean and variances b 2 ff £ SNR respectively and 
Z is statistically independent of {Xi}f =1 with zero mean and 
unit variance. Hence, the minimizing value of 8 is the linear 
MMSE estimation coefficient of X from Y. A straightforward 
calculation shows that the minimizing value of 8 is 

E(XY) _ SNR g T Ba 



8 = 

Var(Y) 

and the MSE it achieves is 



1 + SNR g T Bg 



offf (g. a, B) = mm <7 e 2 ff (g, a,8,B) 



= SNR a T B 



SNR Bgg T B 
1 + SNR ■ g T Bg 



a T (SNR B 



-lTJ-1 



(SNR 



*B 1 



gg ) a 

TN-V2 

gg ) a 



(37) 
(38) 



where again (l38l can be verified using Woodbury's matrix 
identity ||48l Thm 18.2.8]. Accordingly, we define 



^corap(g, a, B) 



1 / SNR 

2 ° g Ue 2 ff(g^,B) 



(39) 



As in Section |nll we define the set of optimal L coefficient 
vectors for the equivalent channel ( |29l as the L linearly 
independent vectors {ai, . . . ,8ll} that yield the highest com- 



putation rates R { 



comp, 1 



^comp(g, ai,B) > ■ ■ ■ > i?, 



comp,L 



-Rcomp(g, a^, B) (see Definition |7j. The compute-and-forward 



transform of the induced L-user MAC is 



melCt bmd-r, 



mod A \ 



^ [PLY + £*=1 a Ll J2 m eK e b mdr, 
t e ff,l \ 

A | : + Z e ff mod A, 

t e ff,L / 



mod A 



(40) 



where A = [ai • • • a L ] T and Z eff = [z? s i ■ ■ ■ z^ s L ] T . 

The next two theorems are simple extensions of Theorems [3] 
and [5] Their proofs are given in Appendix [B] 

Theorem 7: The sum of optimal computation rates for the 
effective L-user MAC ( 1291 is lower bounded by 



L 



RrnmnJ > _ i og __ _ _ _ i og(j L). 



f=i 



Theorem 8: Consider the effective L-user MAC (|29l , in- 
duced from the if -user MAC (|28T >. characterized by the 
effective channel vector g and the effective weight matrix B. 
For any e > and n large enough there exists a chain of n- 
dimensional nested lattices A C C • • • C Ai forming the 
set of codebooks C\, . . . , Cl with rates i?i, . . . , Rl such that 
for all gel 1 and B, if: 

1) each user m £ K,i encodes its message using the 
codebook Hi or a codebook nested in Ci, 

2) the integer-valued matrix from the compute-and-forward 
transform of the effective MAC (|29l can be pseudo- 
triangularized with the permutation vector ir, and the 
optimal computation rates are i? comp ,i > • • • > i? com p ,l, 

3) all rates Ri, . . . ,Rl satisfy 



R( < R c 



for 



1, 



• L 



(41) 



where ir 1 is the inverse permutation vector of it, 
then all effective lattice points t e ff/ can be decoded with error 
probability smaller than e. 

Corollary 2 (Achievable symmetric rate): Consider the ef- 
fective L-user MAC induced from the LT-user MAC d28l . 
characterized by channel coefficients g and the effective 
weight matrix B. There exists a pair of n-dimensional nested 
lattices A C A x forming the codebook L of rate R such that 
for all gel 1 and B, if 

1) all users encode their messages using L (or codebooks 
nested in £), 

2) The Lth optimal computation rate in the compute-and- 
forward transform of (l29i is R CO mp,L, 

3) R < i?comp,L, 

then, for n large enough, all effective lattice points t e ff,^ can 
be decoded with an arbitrarily small error probability. 

Remark 9: Corollary [2] is easily obtained from Theorem [8] 
However, it can also be established without incorporating the 
compute-and-forward transform machinery. Indeed, if all users 
transmit from the same lattice codebook with rate smaller than 
Rcomp,L, by Theorem |2] each of the L equations with optimal 
coefficient vectors can be decoded (without using successive 
decoding as in the compute-and-forward transform approach). 
Then, the decoded equations can be solved for the effective 
lattice points. 

In Section [VI] we introduce two achievable schemes for 
the AT-user Gaussian interference channel. One of them is a 
simple transmission scheme where all users transmit from the 
same nested lattice code. The result of Corollary |2] suffices 
to establish the rates achieved by this scheme. In the second 
achievable scheme, which mimics the Han-Kobabyshi scheme 



SUBMITTED TO IEEE TRANS. INFO THEORY 



14 



for the two-user interference channel, each user transmits 
a superposition of codewords taken from two nested lattice 
codebooks. In this case Corollary [2] does not suffice and 
Theorem [8] which uses the compute-and-forward transform 
machinery, is needed. 

In Section [VT1 we leverage the above results to lower bound 
the capacity of the symmetric Gaussian A-user interference 
channel. 

V. Symmetric Capacity Upper Bounds 

In this section we give upper bounds on the symmetric 
capacity of the symmetric A-user Gaussian interference chan- 
nel. To We follow the same arguments given in J20| for 
showing that the symmetric capacity of the symmetric K- 
user interference channel is upper bounded by that of the 
symmetric two-user interference channel. Namely, eliminating 
all but two users, say users 1 and 2, the symmetric capacity 
is upper bounded by the results of [5 1. This is simply because 
removing interferers cannot decrease the symmetric rates for 
users 1 and 2. Thus, the upper bounds from [5] hold for the 
symmetric rates of user 1 and 2 in the A-user symmetric 
interference channel. Repeating the same argument for each 
pair of users we see that the upper bounds on Csym developed 
in fl5j] for K = 2 continue to hold for all K > 2 as well. 
Therefore, the symmetric capacity of the symmetric A-user 
Gaussian interference channel is upper bounded as |]5] 



C 



SYM 



< < 



'ilog(l 

|log(l- 

|log(l- 
Ulog(l- 



- INR4 
SNR) 
SNR 4 
SNR) 



SNR 
1 + INR 

- 5 log 

INR) 



SNR 
1 + INR 



o< a <§ 

§ < a < 1 

1 < a < 2 

2 < a. 

(42) 



Since we are only after an approximate capacity characteriza- 
tion, we further upper bound Csym for SNR > 1 as 



Cs\ 



< < 



1 + INR 

2INR) 
SNR) 

Thus, for all values of SNR we have 



5 log (2 
±log(l 
|log(l 



2SNR 
1 + INR 

2INR) 

4SNR 2 



< a < \ 



iiog(i 

±log(l 



Csym < < 



\ log 



1 + SNR 
1 ' 1 + INR 

(INR) + 1 
SNR 



1 



INR / 

±log + (INR) + l 
\ log (1 + SNR) 



§ < a < 1 

1 < a < 2 

2 < a. 



< a < \ 



!<«<! 

1 < a < 2 

2 < a. 



(43) 



VI. Achievable Schemes 

This section introduces two simple achievable schemes for 
reliable communication over the symmetric if -user interfer- 
ence channel which are based on nested lattice codes. These 
schemes are then shown to approximately achieve Csym, the 



symmetric capacity of the channel, for all channel gains g, 
except for an outage set of bounded measure. This outage set 
is explicitly characterized. 

We begin by describing the two schemes and deriving their 
achievable symmetric rates. These rates are given in terms 
of the optimal computation rates corresponding to a certain 
effective multiple access channel, i.e., the rates are given 
as a solution to an optimization problem. This optimization 
problem, which amounts to finding the optimal coefficient 
vectors, can be efficiently solved numerically, as described 
in Section IIII-DI Figure [4] shows our achievable rates for the 
three-user symmetric interference channel as a function of the 
interference level g, for several values of SNR. It is evident 
that the obtained rates significantly improve over time-sharing 
even for moderate values of SNR. 

In order to establish the approximate optimality of these 
schemes, we derive explicit lower bounds on the rates they 
achieve which depend only on the SNR and INR. As in the 
two-user case, the symmetric capacity exhibits a different 
behavior for different regimes of interference strength, char- 
acterized by the parameter a. 

We now present the two achievable schemes. The first 
achieves the approximate symmetric capacity in the noisy, 
strong and very strong interference regimes, while the second 
achieves the approximate symmetric capacity in the weak and 
moderately weak interference regimes. 

First scheme - A single-layer lattice code: A pair of nested 
lattices A C Ai is utilized to construct the codebook C = 
Ai n A of rate -Rsym- All users encode their messages using 
this codebook. Since all interferers arrive at the fcth transmitter 
with the same gain, they will be aligned into one effective 
lattice point. Thus, the A-user MAC seen by the kth receiver 
becomes an effective two-user MAC of the form defined in 
Section ITV-CK see Example [3]) 



yfc = x e ff,fei + .gx e ff : /c2 



•*k, (44) 
tm are the effective 



where x eff ,fci = x fe , x effi/c2 = 

users, 6g ff l — 1, b^ s 2 = K — 1 are the effective users' weights 
and g = [1 g] T is the vector of channel gains. 



Transmit 



Equations Decoded by Receivers 



Xl 




anxi +ai 2 y^xi 




a 2 ixi + a 22 ^ xi 








x 2 




anx 2 + ai2 ^ xi 




021X2 + a 22 2_j x ^ 








X_ftT 




anXK + a i2 X£ 




a2iXif + a 22 x^ 



Fig. 8. Illustration of the single-layer lattice scheme. Each transmitter sends a 
codeword drawn from a common lattice. Each receiver decodes two equations 
of the codewords, which it can then solve for its desired message. 

The next theorem gives an achievable rate region for the 
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A-user interference channel when each receiver jointly de- 
codes both the effective user Xgff^i which carries the desired 
information, and the effective user x e ff ^2 which carries the 
sum of interfering codewords. The theorem relies on decoding 
two independent linear combinations of the effective lattice 
points. See Figure [8] for an illustration. This is in contrast to 
the successive decoding technique used in |[T8l , where first 
the interference is decoded and removed, and only then the 
desired lattice point is decoded. 

Theorem 9: Let i? comp ,i > A comp .2 be the optimal com- 
putation rates for the effective MAC (l44l induced by the 
symmetric A-user interference channel ([T}. Any symmetric 
rate Rsym < A CO mp.2 is achievable for the symmetric A-user 
interference channel (|T). 

Proof: Corollary [2] implies that for any symmetric rate 
Asym < -Rcomp.2 there exists a pair of nested lattices A C Ai 
such that both effective lattice points can be decoded at each 
receiver. Since the first effective user x e ff : ki carries all the 
desired information for the fcth receiver, it follows that any 
Asym < -Rcom P ,2 is achievable. ■ 

The next theorem gives an achievable rate region for the 
A-user interference channel when each receiver decodes only 
its desired codeword, while treating all other interfering code- 
words as noise. This theorem can be trivially proved using 
i.i.d. Gaussian codebooks. Nevertheless, we prove the theorem 
using nested lattice codebooks for completeness. 

Theorem 10: Any symmetric rate satisfying 



A SY m < ^ lo § 



1 



SNR 



1 + (A — l)g 2 SNR 



is achievable for the symmetric A-user interference chan- 
nel (D. 

Proof: Decoding x& at the fcth receiver of the sym- 
metric A-user interference channel (Q3, while treating all 
other users as noise, is equivalent to decoding the equation 
with coefficient vector a = [1 0] T in the effective two- 
user MAC ( f44b . Therefore, any symmetric rate satisfying 
Asym < ^comp(g, [1 0] T ,B) is achievable. The effective 
variance for decoding this equation is found using (|37| | to be 



a 2 ff (g, [1 0] T ,B)=SNR 1 



SNR 



1 + (A - l).g 2 SNR 



which, using d39l ), implies that 

A comp (g, [10] T ,B) = ilog(l 



SNR 



1 + (A - 1)</ 2 SNR 



For the two-user case, it is known that in the weak and 
moderately weak interference regimes each receiver should 
decode only part of the message transmitted by the other 
user [5|. A natural extension of this Han-Kobayashi lfl5l 
approach to the A-user case is for each receiver to decode 
linear combinations that only include parts of the interfering 
messages. This is enabled by using a superposition of two 
lattice codewords at each transmitter, as we describe next. See 
Figure [9] for an illustration. 



Second scheme - Lattice Han-Kobayashi: This scheme em- 
ploys a chain of nested lattices A C A2 C Ai to construct 
two codebooks C\ and £2 with rates R\ and i?2, respectively. 
Each user k splits its message Wk into two messages, a public 
message Wki that is mapped into codebook L\ and a private 
message Wki that is mapped into codebook £ 2 . It is convenient 
to treat each user k as two virtual users with codewords x^i 
and Xfe2 that carry messages Wki and Wk2, respectively. User 
k transmits a superposition of its virtual users' codewords, 



Xfc 



Vl-fx H + 7X fe2 , 



for 7 € [0, 1). The signal seen by the fcth receiver is 



Yk 



= V 1 - 7 2 *fci + 7 x fc2 



sVi -7 2 Xml 



91 



x m2 + z fc . 



(45) 



m^k m^k 

which induces the effective four-user MAC 



Yk 



= V 1 - 7 2x eff,fcl + 7 x eff,fc2 



+ 9V 1 - 7 2x eff,fe3 + 57 x eff,fc4 + Zfc, 



(46) 



with effective users x effife i = x fe i, Xeff,fe2 = x fc2, x eff,fe3 = 
E m /t x mi and x eff,fc4 = E m ^ fc x m2 . The effective users' 
weights are 6 2 ff A = 1, 6 2 ff 2 = 1, 6 2 ff 3 = K - 1 and 6 2 ff 4 = 
A - 1, and 



g = y/l - 7 2 7 g\J\ - 7 2 37 

is the vector holding the channel gains. 

The receiver aims to decode the effective codewords x e ff ki, 
x eff.fe2 and Xeff fcs while treating the fourth effective codeword 
x eff.fe4 as noise. The next lemma will be useful for the 
derivation of rates achieved by this scheme. Its proof is given 
in Appendix ICl 

Lemma 1: Consider the effective A-user MAC ( |29| i, where 
the decoder is only interested in the first L—l effective lattice 
points t eff)1 , . . . ,teff,£_i and let k = 1/^1 + SNRg|6 2 
Any rate tuple achievable for decoding t e ff i, . 
the effective (L ~ l)-user MAC 

L-l 

is also achievable for decoding the desired L — l lattice points 
over 1 



,teff,i-i over 



(47) 



The next theorem gives the achievable rate region for the 
lattice Han-Kobayashi scheme. 

Theorem 11: Let k( 7 ) = 1/^1 + SNR ff 2 7 2 (A - 1) and 
consider the effective MAC 



Yk 



= k(7)v /i - 7' 



+ K(7)g\/l - 7"x, n /,,; 



x eff,fcl 

a. 



K(7)7X e ff ife2 
+ Zfc, 



(48) 



with effective channel vector 



«(7)vl ~7 2 «(7)7 K (7)ff\A - T 



and effective users' weights 6 2 ff 1 = 1, fe 2 ff 2 = 1, and 



y eff,3 



A - 1. Let {ai(7),a 2 (7)^3(7)} and A comPi i(7) > 
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Fig. 9. Illustration of the lattice Han-Kobayashi scheme. Each transmitter sends a public (blue) and a private (red) lattice codeword. Each receiver decodes 
three equations of the public codewords as well as its desired private codeword while treating the other private codewords as noise. From these equations, 
the receivers can infer their desired public and private messages. 



-Rcomp,2(7) > -Rcomp,3(7) be the optimal coefficient vectors 
and computation rates. Any symmetric rate satisfying 

-RsYM < max i? CO mp,2(7) + R comp, 6 (7) 
7£[0,1) 

is achievable for the symmetric AT-user interference chan- 
nel CQ. 

Proof: The receiver is only interested in the effective 
lattice points teski, t e ff ) fe2. Nevertheless, we require that it 
decodes the three effective lattice points teff fei, t e ff.fc2 and 
tefffe3- Due to Lemma Q] any rate tuple that is achievable 
over the effective channel d48l > is also achievable for de- 
coding t e fffei, t e ff fc 2 and t e ff fe 3 from the original effective 
channel (f46b induced by the lattice Han-Kobayashi scheme. 

Note that teff^i and t e ff,/c3 are points from the same 
codebook C\ with rate Ry, and t e £f fc2 is a codeword from 
£2 with rate i?2- 

Consider a compute-and-forward transform coefficient ma- 
trix A(7) = [ai(7) &2{l) a 3(7)] T for (EH}. For any such 
full-rank, integer-valued matrix, there exists at least one order 
of pseudo-triangularization. Therefore, there exists a pseudo- 
triangularization of A(7) with at least one permutation vector 

7T. 

Consider first the case where 7r(3) = 2, i.e., the effective 
lattice point t e ff,2 is the last to be removed in the successive 
cancelation decoding procedure of the compute-and-forward 
transform. According to Theorem [H] for any R < i? C omp,2(7) 
and i?2 < -Rcomp,3(7) there exists a chain A C A2 C Ai such 
that t e ff kit t e ff,fc2 and t & s t k3 can be decoded from the effective 
channel (|48T > via the compute-and-forward transform. 

Otherwise, tt(1) = 2 or 7r(2) = 2, which means that the 
effective lattice point t e ff.2 is either removed first or second 
from the proceeding equations in the successive cancelation 
decoding procedure of the compute-and-forward transform. 
According to Theorem [8] for any R < R CO mp,3(j) and 
i?2 < -Rcomp,2(7) there exists a chain A C A2 C Aj such 
that t e ff fci, t e ff ^.2 and t e ff,fc3 can be decoded from the effective 
channel (|48T > via the compute-and-forward transform. 

Since Rsym = Ri + R2, and 7 can be chosen such as to 
maximize -Rsym, the theorem is proved. ■ 



The problem of optimizing the power allocation 7 between 
the private and public codewords, played a major role in the 
approximation of the two-user interference channel capacity 
0. Here, we follow this approach and choose 7 such that, at 
each unintended receiver, the received power of each private 
codeword is equal to that of the additive noise. Specifically, 
in the sequel we set j 2 = l/(g 2 SNR). While this choice 
of 7 may be sub-optimal, it suffices to develop our capacity 
approximations in closed form. The achievable symmetric rate 
for 7 2 = l/(g 2 SNR) is given in the following corollary to 
Theorem [TT] 

Corollary 3: Assume g 2 SNR > 1 and consider the effective 
MAC 



Yk = 




5 2 SNR- 1 
K ■ <? 2 SNR 

,g 2 SNR- 1 
K- 5 2 SNR 



X e ff,fcl 



1 



K ■ g 2 SNR 



X e ff,fc2 



X e ff,fe3 + z k, 



(49) 



with effective channel vector 



/ 5 2 SNR- 1 



1 



'g 2 SNR - 1 



K ■ 5 2 SNR V K ■ 5 2 SNR * V K ■ 5 2 SNR 



(50) 



= 1. and 6 eff.3 
>HK > j^HK 



and effective users' weights b 2 n x = 1, b 2 n 2 

1. Lei \a 1 ,a 2 , a 3 j- anu n comp l ± -ft comp 2 ±_ " comp! 3 
be the optimal coefficient vectors and computation rates for 
this effective MAC. Any symmetric rate 



-Rsym < R, 



HK 

comp, 2 



R 



HK 

comp , 3 



is achievable for the symmetric A-user interference chan- 
nel CD. 

Computing the achievable rates given by Theorem [9] and 
Corollary |3]requires finding the optimal computation rates for 
the effective MACs (l44t and j49l , which involves solving an 
integer least-squares optimization problem (see Section ITH-DI) . 
In the remainder of this section, we derive lower bounds on 
these achievable rates that depend only on the values of SNR 
and INR and can therefore be directly compared to the upper 
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bounds ( f43b . To simplify the exposition, we assume g > in 
the sequel, although all results easily follow for g < as well. 



A. Very Strong Interference Regime 

The sum capacity in this regime was characterized exactly 
by Sridharan et al. [18| using a lattice encoding scheme very 
similar to the one used in Theorem [9] The key difference is 
that in [18] each receiver decodes successively: it first decodes 
the sum of interfering codewords and then subtracts it in order 
to get a clean view of the desired signal. Recall that in our 
scheme, each receiver decodes two linear combinations of its 
signal and the interference. 

We now proceed to lower bound the achievable rate of 
Theorem [9] for the case a > 2, i.e., g 2 > SNR. We do this, 
rather than directly using the results of [18 1, for two reasons. 
The first is that the results of [18| only apply for 



g2> (SNR + 1) 2 



SNR 



(51) 



whereas we need the result to be valid for any g 2 > SNR. 
The second is to show that our lattice encoding and decoding 
framework suffices to achieve the approximate capacity in all 
regimes. 

Using the single-layer scheme presented above, the channel 
seen by each receiver is converted to an effective two-user 
MAC Ml . Applying Theorem [7] to this effective MAC, we 
find that the sum of the optimal computation rates is lower 
bounded by 



^comp.l + -Rcomp,2 > ~ log 



l + SNR(l + g 2 (Jf- 1)) 
K - 1 



- 1. 

(52) 



Let i? comp (g, [0 1] T , B) be the computation rate for decoding 
the lattice equation with coefficient vector a = [0 1] T over the 
effective MAC I® with g = [1 g] and B = diag(l, K — 1). 
The effective variance for the equation a = [0 1] T , which is 
calculated using d37l i. is given in (T53l > at the top of the next 
page. Substituting <7g ff (g, [0 1] T ,B) into (f36t gives 



rp 1 /l + SNR(l + 5 2 (^-l)) 

iWg, [0 If, B) = - log (j ,_\ )(1 g + SNR) 



(54) 



For SNR > 1 we can lower bound 

i? comp (g, [0 1] T ,B) > ilog 



g 2 (K~ 1)SNR 
(K — 1)(1 + SNR) 

>\\og{g 2 )-\. (55) 

We would like to find a lower bound on i? C omp,2- Consider 
two cases. If the coefficient vector a = [0 1] T is optimal, i.e., 

#com P ,i = #com P (g, [0 1] T ,B), we have 

ficomp,2 = ^?comp,l + -Rcomp,2 — ficomp(g 7 [0 B) 

>ilog(l + SNR)-l, 



where we have used d52l and d54T i for the last inequality. If 
the coefficient vector a = [0 1] T is not optimal, we must have 



R 



comp,2 > -Rcomp(g, [0 1] T , B) 



>§ lo s(<7 2 )-~, 



where the last inequality follows from (1551 1. Taking the mini- 
mum of the two bounds above we obtain 



/>W~ inin ( ilog(l + SNR)-l,ilog(«? 2 )-i 



~log(l + SNR)-l, 



(56) 



for any SNR > 1 and g 2 > SNR. Since for any SNR < 1 the 
bound (T56l > is negative, the restriction SNR > 1 is redundant. 
Applying Theorem [9] we conclude that for g 2 > SNR any 
symmetric rate satisfying 



i?SYM < ^log(l + SNR)-l 



(57) 



is achievable, and is within 1 bit of the outer bound (l43l l in 
this regime. 

B. Strong Interference Regime 

The strong interference regime corresponds to 1 < a < 2, 
or equivalently 1 < g 2 < SNR. As in the previous subsection, 
we lower bound R, 



comp,2 



in order to obtain a closed-form 



expression for the achievable symmetric rate. In contrast to 
the very strong interference regime, where the lower bound 
on i? C omp,2 is valid for any g 2 > SNR, here we must exclude 
certain channel gains in order to get a constant gap from the 
outer bound (l43l l. That is, the lower bounds we derive for 
the strong interference regime are only valid for a predefined 
subset of the interval g 2 £ [1,SNR). As we increase the 
measure of this subset, our approximation gap worsens. This 
somewhat strange behavior is to be expected from the existing 
literature. The results of J8] and P31 show that for the A'-user 
interference channel the DoF are everywhere discontinuous. 
The notion of DoF corresponds to a w 1. Since the strong 
interference regime contains values of a near 1, we cannot 
expect to achieve rates which are a constant gap from the 
upper bounds of [5| for all values of g. Instead, we show that 
these upper bounds can be approached up to a constant gap for 
all 1 < g 2 < SNR except for some outage set whose measure 
can be controlled at the price of increasing the gap. We will 
see a similar phenomenon when we analyze the moderately 
weak interference regime. 

Substituting g = [1 g] and B = diag(l, K - 1) into d34l ), 
the optimal computation rate for the effective MAC (l44l can 
be written as 



iWi = ilog(SNR)-ilog(a 2 ) 



(58) 



(7 = mm 

(9,01,02 



({P - ai) 2 + {fig - a 2 ) 2 (K - 1))SNR + A , 

where a 2 is the effective noise variance and the minimization 
is over fi € R, and a = [at a 2 ] G 1? \ 0. From d52l combined 
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IS 



^ ff (g,[0 1] T ,B) = SNR- [ 1 ] 



1 

K- 1 



V 



SNR 


' 1 

K - 1 




' 1 

. 9 


[1 9] 


' 1 

K -I 


1 + SNR[ 1 g ] 


1 

K -I 




' 1 " 

. 9 _ 





SNR • 



(X -!)(! + SNR) 
l + SNR+(if-l) 3 2 SNR 



(53) 



with ( 1581 . we have 

1 /l + SNR(l+.g 2 (if-l)) 



if - 1 



comp,l 



-Rcomp,2 > ^ l0 S 

> ilog( 5 2 SNR) - 1 -i? comp ,i 
= I log (, 9 2 ) + i log (a 2 ) -1, 

and combining with Theorem [9] we see that any symmetric 
rate satisfying 

i?SYM < \ log (g 2 ) + \ log (a 2 ) - 1, (59) 

is achievable over the Jf-user interference channel. Thus, in 
order to obtain a lower bound on Csym it suffices to lower 
bound a 2 . The effective noise a 2 can be bounded as 



a'> min [ [((3 — ax) + ((3g — a 2 ) SNR + /3 . (60) 

/3,01,02 \\ / / 

We first hold /3 constant and minimize over ai,a2. If |/3| > 
l/(2g), the optimal choices for the integers ax,a 2 are 



oi = L/?l ) fl 2 = [fig] ■ 



(61) 



If |/3| < l/(2g), rounding the gains will set both ai and a-i to 
zero, which is not allowed. Since g > 1, the optimal choice is 



a\ =0, a 2 = sign(/3). 



(62) 



Now, we are left with the problem of minimizing (1601 over 
p. Instead of performing this minimization problem, we give 
a lower bound on its solution. We do that by splitting the real 
line into three intervals, and lower bound a 2 for all values of 
j3 within each one. Then, we take the minimum over these 
three bounds. 

Interval 1 : < |/3| < l/(2g) 

In this interval it is optimal to set a 2 = sign(/3). Moreover, 
\/3g\ < 1/2, and therefore \j3g — a 2 | > 1/2. Combining this 
with ( f60b gives 

2 SNR 

o 2 g > —■ (63) 

Interval 2 : l/(2g) < \0\ < 1/2 

Here, it is optimal to set ai = \jf] = 0. Substituting ax = 
in ( f60b gives 



a 2 >/? 2 SNR>^>-^SNR 1 / 2 , 



(64) 



where the last inequality follows since g 2 < SNR in the strong 
interference regime. 
Interval 3 : 1/2 < \j3\ 



Since \(3\ > 1/2, we can always decompose (3 into a sum 
of a nonzero integer q and a real number ip E [—1/2, 1/2) 
giving 



a;> min (<^ + q-ai) 2 SNR 

tp,q,ai,a,2 \ 

+ {qg -a 2 + <pg) 2 SNR + (tp + qf 



> min [ip z + (qg-a 2 + <pgy SNR +2- . (65) 

<p,i, a 2 \ \ / 4 J 

The minimization of ((65) with respect to ip (where the 
constraint ip e [—1/2,1/2) is ignored) can be obtained by 
differentiation. The minimizing value of ip is 

9 



* = 1+9 2 
Substituting tp* into (l65T l gives 

1 



(qg-a 2 ). 



o „ > min 



q-Ai \ 1 + g' 



(qg-a 2 ySNR 



(66) 



which, using the fact that g 2 > 1, can be further bounded by 



a 2 > iminmax ( -z(qg - a 2 ) 2 SNR,<f 

" 4 q,a 2 \g z 



(67) 



We would like to obtain a lower bound on a 2 which is valid for 
all g ^ S, where S is an outage set with bounded measure. 
Assume first that g E [6,6+1) for some integer 1 < 6 < 
%/SNR. Let 

g max>6 4 1 SNR 1 / 4 -^ 2 , (68) 
y/b + 1/2 

for some (5 > to be defined later, and note that <7 maX! 6 is not 
necessarily an integer. Define 

$ b = ^6+1/2 SNR- 1/4 " 5/2 (69) 

and let 5b be the set of all values of g E [6. 6 + 1) such that 
the inequality 

\qg-a 2 \<<S> b (70) 

has at least one solution with integers q and a 2 , where q is in 
the range < q < q m ^ b . Let S b = [b,b+l)\S b . By <|67]), <|68l> 
and ( [69) . we have that for all g E S b 

a 2 > — min [ min max [ — (qg — a 2 ) 2 SNR, q 2 ) , 

" 4 V <9<L9n, 1 »,i,J,Q2 V5 " / 



min max — (qg ~ a 2 ) SNR, q' 

\q«^.b\>q,a. 2 \g z 



1 




> - 


min 


- 4 




1 






min 


~ 4 





I 2 

imax,6 



^iZ^SNR^.^^SNR 1 / 2 -^. (71) 

g 2 '6+1/2 ; 
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%(2) |1 



6 + 1 



-q 



for all y S [6,6+1) up to an outage set Sb of measure not 
greater than 2~ c . 

Combining the three bounds (|63k (l64l and (|76T > gives 



6+1 



r 6 (3) LL 



□ 



n 



6+i 



6+| 



-q 

6+1 



5„ 



P- 



n □ n 



i r 



6+i 6+1 6+§ 



6+1 



Fig. 10. An illustration of the sets 7{,(1), 7f,(2), 7{,(3) and their union St- 
In this illustration g m!1Xi i, = 3 and <E>f, = 1/16. 



Since 6 > 1, we have that 

§<6+i<2 ff . 
Thus, d76l > can be further bounded by 

a 2 n > — ^SNR 172 - 5 . (72) 

We now turn to upper bound the Lebesgue measure of the set 
Sb- Let X = [—1,1) and define the set 



76(g) 



6,6+-, 
q 



9 



—I 

q . 



mod [6,6+1), 



where the sum of the two sets is a Minkowski sum. Writ 
ing ([70]) as 



a.2 

.9 

q 



< 



q 



(73) 



we see that for a given q and g G [6,6+1) the inequality 
admits a solution if and only if g G %{q)- It follows that 



S b = U T "(«)- 

9=1 



(74) 



See Figure [H)]for an illustration of the sets %{q) and Sb- Thus, 
the Lebesgue measure of Sb can be upper bounded by 

fi{S b ) = Vol {S b ) 

Lfc,x. 6 J 
< ^ Vol(T 6 (g)) 



9=1 

La™*,*] 



< 



9=1 

= 2SNR 



-<5 



(75) 



Setting <5 = (c + l)/log(SNR) and substituting into (1721 
and d75l ) gives 



fj? > ^SNR 1 / 2 
9 " 16 fll 



(76) 



CT 2 > min f SNR ( ^SNR 1 / 2 , ^SNR 



1/2 



> 



16| 5 



V 4 '4| 5 | 
■SNR 1/2 



16 «7 



(77) 



for all g £ [6, 6 + 1) up to an outage set Sb of measure not 
greater than 2~ c . 

Combining d59l > and ( TTTb we see that for all g G [6, 6 + 1) 
up to an outage set Sb of measure not greater than 2~ c any 
symmetric rate satisfying 



i?SYM<^log( 5 2 SNR)-^-3 
= ilog(INR)-|-3 



(78) 



is achievable. We conclude that the symmetric rate (fTFt is 
achievable for all channel gains in the strong interference 
regime except for an outage set whose measure is a fraction 
of 2~ c of the interval 1 < \g\ < VSNR, for any c> 0. 

Remark 10: We note that for any c > the set of channel 
coefficients that fall in the outage set S can be easily deter- 
mined. 

C. Moderately Weak Interference Regime 

The moderately weak interference regime is characterized 



by 2/3 < a < 1, or equivalently, SNR 



-1/3 



< g z < 1. As 



in the strong interference regime, we show the achievability 
of symmetric rates which are a constant gap from the upper 
bound for a certain fraction of the channel gains. As opposed 
to the very strong and strong interference regimes, where a 
single-layered lattice scheme suffices to achieve the approxi- 
mate capacity, here we will need the second scheme, which 
employs two layers of lattice codes at each transmitter. 

As mentioned in the description of the second scheme, 
we set the power of the private lattice codewords so that 
they are perceived at noise level at the unintended receivers. 
The achievable rate for this choice is given by Corollary [3] 
and we now turn to lower bounding this achievable rate in 
closed form. Let g\, <?2 and 53 be the channel gains in the 
effective three-user MAC (|49T > from Corollary [3] and recall 
that, for this effective channel, the effective weight matrix is 
B = diag(l, 1, K — 1). We begin by applying Theorem [7] to 
the effective channel d49b to get the following lower bound on 
the sum of optimal computation rates, 



HK 

comp,z 



1 ( 1 + SN R(g\ + g\ + (K - 1) 9 |) \ 8 
~ 2 ' 



K - 1 



; 1 -^og(3) 



> 1 log (SNR( 3 2 + gl)) - 1 log (3 3 (K - 1)) 
= l - log (SNR)-i log {27 K(K- I)). 



(79) 
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The optimal computation rate can be written as 



R, 



1 



1 



log(SNR)--log(a 2 K ) 



comp,l 2 v~" "* v 2 



(80) 



where er^ K is given in ( f8TT > at the top of the next page. The 
minimization in (fSTb is performed over all j3 € K and a™ = 
[a\ a,2 a 3 ] e Z 3 \ 0. Applying Corollary [3] we see that any 
symmetric rate satisfying 

Rsym < \ log(^HK) - \ log {17K{K - 1)) (82) 

is achievable for the if -user interference channel. Therefore, it 
suffices to lower bound the effective variance er^ K . Substituting 



Now, the problem of lower bounding <t^ k reduces to min- 
imizing (f84t over f3. Rather than solving this cumbersome 
minimization problem, we split the real line into four intervals, 
and lower bound cr^ K for all values of f3 within each one. 
Then, we take the minimum over these four lower bounds. 
In a similar manner to the previous subsection, we define 
5 = (2c + 8)/ log(SNR), where c > is some constant. The 
lower bounds below are derived in Appendix ID-AI 



Interval 1 : < 



< 1/2 



Kg*Sm - 
/3 -V<7 2 SNR-l /3 

in f8ll , which is allowed since f3 can take any value in K, 
gives 



^HK > 



SNR 



(87) 



Interval 2 : 1/2 < |/3| < v^fSNR 1/4 - a/2 /2 



For all values of SNR < \g\ < 1 except for an outage 
set with measure not greater than 2~ c we have 



mm 

/3,oi,o 2 ,o 3 



- 2 Kg 2 SNR 
' 5 2 SNR- 1 



+ (/3-ai) 2 SNR 



> 



2- 2c SNR 1/2 



4 -2 s 



(88) 



-^=====-a 2 \ SNR+(^-l)(/3 5 -a 3 ) 2 SNR). interval 3 : ^[SHR 1 '^' 2 /2 < |/3| < SNR 1 / 4 / -\/8|g| 



(83) 

In the sequel, we assum^ SNR > 4. With this assumption, 
Vg 2 SNR- 1 > 1 for all g 2 > SNR~ 1/3 , i.e., for all values 
of g in the moderately weak interference regime. We will also 
use the fact that the inequality ^/g 2 SNR — 1 > 1 continues 
to hold for all g 2 > SNR -1 / 2 , i.e., for all values of g in the 
weak interference regime. Since ^/g 2 SNR - 1 > 1 implies 
5 2 SNR > 2, we also have 3 2 SNR/( 3 2 SNR - 1) > 1 and 
hence d83l can be lower bounded as 



^HK > 



Interval 4 : SNR 1/4 /^8|y| < 



SNR 



1/2 



4 -2 s 



(89) 



^-2 \ 



1 SNR 



1/2 



(90) 



'HK 



> _ min K0 2 + 0- ai) 2 SNR 



/3, oi,o 2 ,o 3 



^ 5 2 SNR- 1 



Combining the four lower bounds d87), (|88), d89]l and d90j> 

we have 

? /l_ 2~ 2c SNR 1/2 1SNR 1/2 \ 

a HK > min I -SNR, 



a 2 j SNR+(^-l)(/3 ff -a 3 ) 2 SNRj. 

(84) 



.4" ' 2 10 
2- 2c SNR 1/2 



3 2 '8 



We first hold f3 constant, and minimize over 01, ci2, (X3. If > 
1/2, the optimal choices for the integers a\, a 2 , a 3 are 



! = [ft, a 2 = ^/V.g 2 SNR-ll , a 3 = [fig] 



(85) 



If |/3 1 < 1/2 all three integers ai,a2,a 3 from (185l l are zero 
and thus are invalid. Therefore, for these values of j3 one of 
the integers must take the value 1 or —1. Since for SNR > 4 
and SNR~ 1/2 < g 2 < 1 we have 

max |/3/vVSNR-l|, \p g \) = 



for all SNR -1 / 3 < g 2 < 1 up to an outage set of measure not 
greater than 2~ c . Thus, substituting our lower bound for <t hk 
into (1821 . we find that any symmetric rate satisfying 



1 / cm pi/ 2 \ 1 1 

i?SYM < 7T log f - c - 5 - - log(27) - - log(if 2 ) 



the optimal choices of ai,a2,a 3 for values of |/3| < 1/2 are 
ai = sign(/3), a 2 = 0, a 3 = . (86) 

6 This assumption is valid, since for SNR < 4 the symmetric capacity is 
upper bounded by l/21og(l + 4) = 1.161bits. Our capacity approximations 
in this subsection, and also in the next subsection, exhibit a constant gap 
greater than 7/2 bits, and therefore hold for SNR < 4. 



<\ log 



'9 

'SNR 1/2 



'9 
SNR 

x/TNR 



j-c-8- log(iv-) 
- c - 8 - log(if) 



(91) 



is achievable over the symmetric if-user interference channel 
for all SNR -1 / 3 < g 2 < 1 up to an outage set of measure not 
greater than 2~ c . 

7 If ^/[y|SNR 1 / 4 " 4 / 2 /2 < 1/2 this interval is empty, and we skip to 
interval 3. 
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(81) 



D. Wea£ Interference Regime 

This regime is characterized by 1/2 < a < 2/3, or equiva- 
lently, SNR~ 1/2 < g 2 < SNR~ 1/3 . As in the moderately weak 
interference regime, we develop a closed-form lower bound on 
the achievable symmetric rate of Corollary A key difference 
is that the bound derived here is valid for all channel gains, 
rather than up to an outage set. 

We first note that equations (|82l and (|84l continue to hold 
in this regime as well as in the moderately weak interference 
regime, and the optimal choices of 0,1,0,2,0,3 are also as in (l85T l 
and (I86t , As before, we divide the real line into four intervals, 
give lower bounds on cr HK which hold for all values of j3 in 
each one, and conclude that er HK is lower bounded by the 
minimum of these four bounds. The lower bounds below are 
derived in Appendix ID-BI 

Interval 1 : < |/3| < 1/2 



^HK > 



SNR 



Interval 2 : 1/2 < |/3| < l/{2\g\) 



CT HK > 



5 2 SNR 



Interval 3 : 1/(2| 5 |) < < ^.g 2 SNR/8 



ct hk > 



1 



(92) 



(93) 



(94) 



Interval 4 : ^ ff 2 SNR/8 < 



2 .g SNR 
ct hk > — ; — ■ 



(95) 



Combining the four lower bounds (j92), (J93j, dHJ and (|95t 

we have 



~2 ^ 



1 



4 

ff 2 SNR 



min (SNR, 3 2 SNR, 5 - 4 



(96) 



where d%} is true since SNR~ 1/2 < g 2 < SNR _1/3 . It follows 
by substituting d96*i l into (l82l that any symmetric rate satisfying 



1. fg 2 SNR 

i?SYM < g l0 § I 4 



ilog(27)-ilog(X 2 ) 



<ilog ( 5 2 SNR)-^-log(K) 



Ilog(INR)-^ 



log(X) 



(97) 



is achievable for the symmetric K-nsei interference channel 
with SNR~ 1/2 < q 2 < SNR" 1/3 . 



E. Noisy Interference Regime 

The noisy interference regime is characterized by a < 1/2, 
or equivalently g 2 < SNR -1 / 2 . In this regime, each receiver 
decodes its desired codeword while treating all interfering 
codewords as noise. Lattice codes are not necessary in this 
regime in order to approximate the symmetric capacity: ran- 
dom i.i.d. Gaussian codebooks suffice. Nevertheless, the same 
performance can be achieved with lattice codes as shown in 
Theorem [10] which states that any symmetric rate 



-RSYM < ^ log ( I 



SNR 



1 + (K — 1)<7 2 SNR, 
is achievable. It follows that any symmetric rate satisfying 

1 



1. A SNR 

i?SYM < 2 log I 1 + - 



1 



log 1 



+ .g 2 SNR 
SNR 



1 + INR 



1 



^log(iC-l) 
\og(K - 1) 



(98) 



is achievable. 



VII. Degrees-of-Freedom 



We have shown that the compute-and-forward transform can 
approximate the capacity of the symmetric X-user interference 
channel up to a constant gap for all channel gains outside 
a small outage set. Ideally, we would like to use a similar 
approach to approximate the capacity of the general (non- 
symmetric) interference channel. In contrast to the symmetric 
case, where all interferers are automatically aligned (if they 
all use the same lattice codebook), the interferers will be 
observed through different channel gains. A linear combination 
of lattice codewords is a codeword only if all of the coefficients 
are integers. Thus, in order to induce alignment, all of the 
interfering gains should be steered towards integers, which is 
an overconstrained problem. 

Recently, Motahari et al. [ 8 ] introduced a precoding scheme 
which achieves K/2 DoF for almost every channel realization 
of a JC-user Gaussian interference channel. Their technique 
relies on symbol-level alignment. Namely, each user encodes 
its message into several layers. Each symbol in a layer belongs 
to the one-dimensional lattice Z. Each user superimposes its 
layers in a clever way, such that at each unintended receiver 
symbols from different layers are received aligned with those 
transmitted by other interfering users. The resulting effective 
channel at each receiver is an effective multiple access channel 
of the type introduced in Section IIV-CI When the numbers 
of layers is large, roughly half of the effective users at each 
receiver carry valuable information while the other half are 
effective interferers. Through the use of techniques drawn 
from Diophantine approximation, they showed that each of 
the L effective users achieves essentially 1 /L DoF. Thus, each 
receiver obtains 1/2 useful DoF. 
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In this section, we show that the same asymptotic results 
can be achieved using the compute-and-forward transform. 
Specifically, we show that nested lattice codes allows each 
of the L effective users in the effective MAC ( f29b introduced 
in Section IIV-CI to achieve X/L degree of freedom for almost 
all channel realizations. Therefore, the one-dimensional lattice 
used in HI can be replaced with 71-dimensional lattice codes 
while maintaining the same asymptotic performance. The ben- 
efit of using nested lattice codes combined with the precoding 
scheme of [8| is that one can obtain an achievable rate region 
at any SNR rather than only a DoF characterization. We begin 
by analyzing the DoF offered by the if optimal computation 
rates in the if -user MAC. Then, we translate the results to the 
effective L-user MAC. 

Theorem [3] in Section |IV] guarantees that the sum of the 
optimal computation rates is close to the sum capacity of the 
MAC. However, the theorem does not tell us how the sum rate 
is divided between the if rates. We now show that, in a DoF 
sense, the sum is equally split between all if rates for almost 
every channel realization. Recall the definition for DoF: 

-Rcomp,fc(SNR) 



J'Comp.A; 



= lim 

SNR->oo i log(l + SNR) 



(99) 



First, we upper bound d C omp,i, the DoF provided by the best 
lattice equation. 

Theorem 12: Let fx, ■ ■ ■ , fx be functions from R m to R 
satisfying 

1) f k for k = 1, . . . , K is analytic in R m , 

2) 1, fx, . . . , fx are linearly independent over R, 
and define the manifold 



M = {[fx(h) ■■■ f K (h)] : hel m } 



(100) 



For almost every h E Ai the DoF offered by the best lattice 
equation is upper bounded by 

1 



^comp.l 



< 



if 



(101) 



The proof is given in Appendix [E] and is based on showing 
that restricting the scaling coefficient j3 from (0 to the form 
j3 = q/hx for q € Z incurs no loss from a DoF point of view. 
This way, the first coefficient of j3h is an integer. Then, a result 
from the field of Diophantine approximation which is due to 
Kleinbock and Margulis | 57] is used in order to lower bound 
the error in approximating the remaining K — 1 channel gains 
with integers. 

Remark 11: The manifold Ai from Theorem Q~2] which 
appears in all of our DoF results in this paper, is the same 
manifold used in [8| to prove the achievability of if/ 2 DoF 
for the if -user interference channel. Thus, all achievable DoF 
results from [8| are also achievable using the compute-and- 
forward transform. 

As a special case of Theorem [T2] we may choose the 
manifold Ai as R A which implies the following corollary. 

Corollary 4: For almost every h g R^ the DoF offered by 
the best lattice equation is upper bounded by 

1 



Remark 12: Niesen and Whiting [58 1 studied the DoF of- 
fered by the best lattice equation and showed that 



^comp.l — 



1/2 if = 2 
2/(K + l) if>2 



(102) 



for almost every h g R^. Our bound therefore agrees with 
that of J58) for if — 2 and improves it for if > 2. 

The next corollary shows that all if optimal computation 
rates offer 1/if DoF for almost every h satisfying mild 
conditions. 

Corollary 5: Let Ai be a manifold satisfying the conditions 
of Theorem [12] For almost every h £ Ai the DoF provided 
by each of the if optimal computation rates is d CO mp,k = 1/if. 

Proof: Theorem |3] implies that 2~2k=i ^comp,fc > 1- Using 
the fact that d comPj fe is monotonically decreasing in k and that 
c^comp.i < f° r almost every h e Ai, the corollary follows. 

■ 

The corollary above implies that in the limit of very high 
SNR not only is the sum of computation rates close to the sum 
capacity of the MAC, but each computation rate scales like the 
symmetric capacity of the MAC for almost all channel gains. 
Note that our analysis (as well as that of [58|) is within the 
context of the achievable computation rates stemming from 
Theorem [2] 

The next corollary follows from Corollary |3]and Theorem 

Corollary 6: Let Ai be a manifold satisfying the conditions 
of Theorem Q~2] The DoF attained by each user in the K- 
user MAC under the compute-and-forward transform is 1 /if 
for almost every he Ai. In particular, the DoF attained by 
each user in the if -user MAC under the compute-and-forward 
transform is 1 /if for almost every h £ R*" 

The next theorem shows that for almost every effective 
L-user multiple access channel of the form introduced in 
Section IIV-CI each of the effective users achieves 1/L degree 
of freedom. The proof is given in Appendix [F] 

Theorem 13: Let fx, ■ ■ ■ , Jl be functions from R m to R 
satisfying 

1) f e for I = 1, . . . , L is analytic in R m , 

2) 1, fx, ■ ■ ■ , ft are linearly independent over R, 
and define the manifold 

M = {[fx(g) ■■■ f L (s)] ■■ geR m }. 

For almost every g e Ai the DoF offered by each of the L 
optimal computation rates for the effective MAC (|29| | is 

-Rcomp^(SNR) _ 1 

~ L' 



^comp,£ 



lim 



SNR x> ilog(l + SNR) 



VIII. Discussion 



(103) 



^comp , 1 



< 



if 



In this paper, we have developed a new decoding framework 
for lattice-based interference alignment. We used this frame- 
work as a building block for two lattice-based interference 
alignment schemes for the symmetric real Gaussian if -user 
interference channel. These schemes perform well starting 
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from the moderate SNR regime, and are within a constant gap 
from the upper bounds on the capacity for all channel gains 
outside of some outage set whose measure can be controlled. 

A natural question for future research is how to extend 
the results above to the general Gaussian K-usei interference 
channel. The main problem is that, in the general case, the 
interfering lattice codewords are not naturally aligned, as their 
gains are not integer-valued. Therefore, in order to successfully 
apply lattice interference alignment, some form of precoding, 
aimed towards forcing the cross channel gains to be integers, 
is required. Unfortunately, simple power-backoff strategies do 
not suffice, even in the three-user case. 

One option for overcoming this problem is to use many 
layers at each transmitter, as in [8 1, and create partial alignment 
between interfering layers. While this achieves the optimal 
DoF, it performs poorly at reasonable values of SNR, as there 
will be a rate loss for each additional layer. As a result the 
rate region obtained by combining the compute-and-forward 
transform with the precoding scheme of B is inferior to 
that obtained by time-sharing, for values of SNR of practical 
interest. Another option is to precode not only using power- 
backoff, but also over time, which may partially compensate 
for the lack of sufficient free parameters. An example for such 
a precoding scheme is the power-time code introduced in [21 1. 

A positive feature of the compute-and-forward framework 
is that it does not require perfect alignment of the lattice 
points participating in the equations. Namely, the effect of 
not perfectly equalizing the channel gains to integers is an 
enhanced effective noise. For the general interference channel, 
this suggests that it may suffice to find precoding schemes that 
only approximately force the cross-channel gains to integers. 

Another issue left for future work is extending the ca- 
pacity approximation for the real symmetric Gaussian K- 
user interference channel to the complex symmetric case. 
The achievable rate region, as given by Theorems [9] and 
Corollary [3] naturally extends to the complex case using the 
results of [19] for the complex case. An open question is 
whether the lower bounds on these achievable rate regions, 
which depend only on SNR and INR, also extend to the 
complex case. The answer seems positive since the bounds 
mainly rely on volume considerations that seem to remain 
valid in the complex case. 

Appendix A 
Proof of Theorem[5] 

We begin with two lemmas which will be useful for the 
proof of Theorem 

Lemma 2: Let A be a K x K matrix with integer entries 
of magnitudes bounded from above by some constant a max . If 
there exists a lower triangular matrix L with unit diagonal such 
that A = LA is upper triangular up to column permutation 
7r, then for a prime p > K (Kl) 2 (Ka max ) 2K a max there also 
exists a lower triangular matrix L' p ) with elements from 
{0, l,...,p-l} and unit diagonal such that AW = [L^A] 
mod p is upper triangular up to column permutation it. 

Proof: Assume that there exists a lower triangular matrix 
L with unit diagonal such that A = LA is upper triangular 



up to column permutation w. We begin by showing that all 
elements in the ith (i > 1) row of L can be written as rational 
numbers with the same denominator 1 < < Kl(Ka 2 a . dx ) K . 
To see this note that if A is triangular up to column permu- 
tation vector 7r, then its ith row contains at least i — 1 zeros, 
namely ay = for j = ir(l), . . . , tt(% — 1). Since L is lower 
triangular, the following equations must hold 

i 

aij = ^2 £ima m3 = 0, for j = tt(1), ...,n(i- 1). (104) 

m— 1 

By definition in — 1, therefore (11041 i can be written as 

i-l 

ii m a m j = -ay, for j = 7r(l), . . . , n(i - 1). (105) 

m— 1 

Define the vectors i^ = [in ■ ■ ■ i^i-\] T , 

a^, 71 ") = _[ a . .j. . . . a i n(i-i)] T and the matrix 

(ai7r(i) ■•■ a i-i7r(i) \ 
a lir(i-l) ■ • ■ ai-lir(i-l) / 

We have, 

A (i,^W = a (v0. (106) 

From the fact that A can be pseudo-triangularized with per- 
mutation vector 7r we know that the system of equations (1106b 
has a solution. Assume that 

rank (a^) = u < K. 

It follows that there are u linearly independent columns in 
a (i,tt) 14 G {1, ... , K} be a set of indices corresponding 
to u such linearly independent columns, and U be its comple- 
ment. Let Ay € Z Axu be the matrix obtained by taking 
the columns of A^ 1 ' 77 ' with indices in U. Since (1106b has a 
solution, we have a^ ,7r ' g span ^A^' 71 ^ . Thus, we can set 

(k) = for all k G U, and (1106b will still have a solution. 
Letting t$ 6 M ,lXl be the vector obtained by taking from £^ 
only the entries with indices in U, it follows that 

A$ r) $ - a(^> (107) 

has a solution. Now, multiplying both sides of ( 1107b by 

( A u* ] ) 8 ives 

A' (M °4 l) =a'< i ' w \ (108) 

where A' ih7r) = (a^' 70 ) T A ( ^ ] £ Z uxu is a full-rank 

matrix and a'^ 7 ^ = (A^ ,7r) ) T sS^ € Z" xl . Note that 

all entries of A.'^ 1 ' 77 ' as well as all entries of a'^' 7 ^ have 
magnitude bounded from above by a max = ua^ ax . We have 

Letting qi = | det(A'^ :7r ')| and applying Cramer's rule for 
matrix inversion (see e.g. fl"8l ) we see that all elements of 
iff can be expressed as rational numbers with denominator qi. 
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Recall the Leibnitz formula (see e.g., [48 1) for the determinant 
of annxn matrix G 



det(G)= sign^Il G ' 

aes. 



(109) 



where S n is the set of all permutations of {l,...,n}. It 
follows that det (A'( l,7r )) must be an integer and in addition 
1 < | det (A'fc*)) | < u\{a mm y. Thus, 1 < Qi < u\(h m . dX ) u . 
Moreover, Cramer's rule also implies that the numerator of 
each element in 6jj is an integer not greater than u!(a max )" 
in magnitude. Since u < K, and since each element of 
is either zero or corresponds to an element in iff, 
each element lij, j < i of L can be written as a rational 
number £jj = rriij/qi with 1 < < if^ifa;^) and 
\m ij \ < K\{Kal !& ) K for » _= 1,...,A\ 

Now, define the matrix L = diag(qi, . . . , q_ff)L and note 
that L e Z KxK due to the above. Let A' (p) = [LA] mod p. 
Since multiplying a row in a matrix by a constant leaves its 
zero entries unchanged, the entries of the matrix 

A'(p) 



[LA] mod p 
= [diag(gi, . . . , g*r)LA] mod p 
= [diag(gi , • • • , qx) A] mod p, 



(HO) 



are zero whenever the entries of A are equal to zero. More- 
over, since all elements of L are bounded in magnitude by 
Kl(Ka^ nax ) K and all elements of A are bounded in magnitude 
by a mdx , all elements of A = LA are bounded in magnitude 



by KK\(Ka\ 



\X, 



qi <K\{Kal^ K , we have 



Combining with the fact that 1 < 



*max 

for all i = 1, . 



< K(K\) 2 {Ka n 



\2X, 



. , K, j = 1,,.,,K. Therefore, for a prime 



number p > K(Kl) 2 (Ka m:lx ) 2K a max the modulo reduction 
in ( II 10b does not change any of the non-zero entries of 
diag(qi, . . . , qx)A to zero. 

Recall that if A can be pseudo-triangularized with a matrix 
L and permutation vector tt then <Zj )7r (i) ^ 0, and hence 
also a.,.s ^ for i = 1,...,K. We have therefore 
shown that for p large enough there exists a lower-triangular 
matrix LAW = [diag(gi, . . . , 5ic)L] mod p with elements from 
{0, 1, ... ,p - 1} such that A'W = [Lb) A] mod p is upper- 
triangular up to column permutation tt. In order to complete 
the proof, it is left to transform L' p ) to a lower-triangular 
matrix with elements from {0, 1, . . . ,p— 1} and unit diagonal. 
Let (<7i) -1 be an integer that satisfies [(qi)~ 1 qi\ mod p = 1. 
Such an integer always exists since is an integer different 
than zero, and p is prime. It is easy to verify that the matrix 
Jj(p) = [diag ((qi)" 1 , ■ ■ ■ , {qx L^] mod p is a lower- 
triangular matrix with elements from {0,l,...,f> — 1} and unit 
diagonal, and A^ p ' = [L^A] mod p is upper triangular up 
to column permutation tt. ■ 

Lemma 3: Let ti,...,tfc be lattice points from a chain 
of nested lattices satisfying the conditions of Theorem 
Let v = EfcLi fl fctfc] mod A and u = Et=i mod A 
be integer-valued lattice equations of these points. Then 

" K 

^ ((«fe + frfc) modp)t fc 
.fc=i 



[v + u] mod A 



mod A. 



Proof: Due to the distributive property of the modulo 
operation we have 



[v + u] mod A 



K 

E 

k=l 
K 

E 

Lfe=l 
K 

E 

if 

fe=i 



a k + b k )t k 



mod A. 



[(fflfc + &fc) mod p + M fe • p]t fc 



((afc + 6fc) mod p)t* 



mod A 



[p • tfe] mod A 



mod A 



(HI) 



where {M k } K 



[P ■ t* 



fe= 
mod A 



are some integers. Utilizing the fact that 
for all lattice points in the chain, which 
follows from Theorem l2]|c) . the lemma is established. ■ 

We are now ready to prove Theorem 



Proof of Theorem \5$ Let T = [tf 



V = 



and 

vjj ] T = [AT] mod A. The compute-and- 



forward transform of the MAC (0 can be written as 







( tl ^ 






s = 






+ z eff 


mod A 













[AT 

[V + 



f Z e ff] mod A 
Z e ff] mod A. 



Assume there exists a pseudo-triangularization of A with 
permutation vector tt, i.e., there exists a lower triangular matrix 
L with unit diagonal such that A = LA is upper triangular up 
to column permutation tt. Lemma [2] implies that there exists a 
lower triangular matrix L( p ) with elements from {0,1, ... ,p— 
1} and unit diagonal such that = [L^ p 'A] mod p is 
upper triangular up to column permutation tt. Since L^ p ^ has 
a unit diagonal it can be written as L'^ = I + R where I is 
the identity matrix and R has non-zero entries only below the 
main diagonal. 

Assume the receiver had access to the side information 
Vi, . . . , \k-\- As the entries of R are non-zero only below 
the main diagonal, the receiver can compute R • V, add it to 
S and reduce modulo A, giving rise to 

S SI = [S + R • V] mod A 

= [AT + RAT + Z ef f] mod A 
= [(I + R)AT + Zeff] mod A 

mod A 



L (p) AT + Z ef f 
[L (p) A] mod p 



Jeff 



mod A 



(112) 





f U \ 






A Cp) 




+ Z e ff 


mod A 











where (TTT21 follows from Lemma[3] Let V = [A^T] mod A 
and recall that is upper-triangular up to column permuta- 
tion tt, thus Oj^j.) = for all j = 7r(fc) + 1, . . . , K. It follows 
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that for any k < K the lattice point t^) does not participate 
in any of the lattice equations Vfc+i, . . . , vr- 

Assume the mapping function between users and lattices 
is chosen as 0(k) — Tr~ 1 (k), i.e., each user k employs the 
codebook Ck = A^-ia) n A. In this case, the densest lattice 
participating in lattice equation is A/.. The decoder uses 
sf 1 in order to produce an estimate 

v fc = [Q Afc (sf)] mod A 



(113) 



for each one of the lattice equations It follows from 
Theorem [2] that there exists a chain of nested lattices A C 
A-k C • • • C Ai forming the set of codebooks Ci, . . . ,Ck 
with rates i?i , . . . , Rk such that all equations vj. , . . . , vjf can 
be decoded with a vanishing error probability as long as the 
rates of all users satisfy the constraints of ((25}. 

We have shown that if the receiver has access to 
Vi, . . . , vk-i it can decode the set of equations V. We now 
show a sequential decoding procedure which guarantees that 
the receiver has the right amount of side information at each 
step. First, note that 



-SI 



fc-1 

E 

m— 1 



mod A, 



(114) 



thus the necessary side information for decoding is only 
Vi, . . . , Vfc_i. In particular, s^ 1 = si and hence vi can 
be decoded with a vanishing error probability with no side 
information. After decoding vi the receiver has it as side 
information, and can therefore compute s^ 1 and decode V2. 
As V2 = [f2iVi + V2] mod A and the receiver knows Vi, it 
can use it in order to recover V2. Now, the receiver has vi 
and V2 as side information and can use it to compute S3 1 . The 
process continues sequentially until all equations vi , . . . , Vk 
are decoded. 

When the process ends we are left with K noiseless 
equations 



vi 





f *1 1 




aw 


t 2 











mod A. 



(115) 



Since A' p ) is upper-triangular up to column permutation, and 
in particular full-rank modulo p, the original lattice points 
ti, . . . , tx each user transmitted can be recovered. ■ 

Appendix B 
Proof of Theorems[7Jand[8] 

Proof of Theorem 1 The proof is identical to that of 
Theorem with F = (SNR _1 B _1 + gg T )~ 1/2 . ■ 

Proof of Theorem |S} Let 







( t e ff4 > 






s = 


} 




+ z eff 


mod A 






\ t e ff,L / 







be the compute-and-forward transform of the effective L-user 
MAC, and assume that A can be pseudo-triangularized with 



permutation vector tt. Repeating the proof of Theorem it 
easy is to see that for any set of rates 



Re < R c 



1,. 



L 



, 'Comp,7r _1 (I) ) 

there exists a chain of nested lattices AC A^ C • • • C Ai 
inducing the codebooks Cg = A^-im HV with rates Re, such 
that if t e ff.^ € Ci for all £ = 1, all effective lattice 

points can be decoded from S. 

If each of the users m G K,g that comprise effective user £ 
uses the lattice codebook Cg (or any codebook nested in Ci), 
then t e ff,^ G Ci and all effective lattice points can be decoded. 



Appendix C 
Proof of Lemma[T] 

In order to decode the desired effective lattice points, it 
suffices to decode L — 1 linearly independent lattice equations 
of them, in which t e ff,i does not participate. Let a = 
\fli ■ ■ ■ a_L_i 0] be some coefficient vector for such an 
equation. The effective rate for computing the lattice equation 
v = [J2e=i ^t^eB,e] mod A with the coefficient vector a over 
the channel d29l ) is 

1, ( SNR 



R, 



comp 



;g,a,B) 



2 l0g V^ 2 ff (g ; a,B) 



(116) 



where 



L-l 



a 2 eff (g, a, B) = min SN R V 0g e 

+ ^(H-SNR ff !6 e \ L ) 

L-l 



= minSNRVV/3/t^ - a e ) 2 b t 



2 h 2 
eff,( 



+ /3 2 , (117) 



where (II 171 i follows by substituting j3 = [3k. The effective 
variance and computation rate for computing an equation with 
coefficient vector a = [ai ••• ci£-i 0] over the effective 
channel d29b are therefore the same as those of computing 
an equation with a = [ai ■ ■ ■ over the effective 

channel d47b . Thus, for purposes of computing equations with 
the effective lattice points t e ff,i, . . . , t e ff,i-i the two channels 
are equivalent. Since this is all we need in order to decode 
t e fr,i, . . . ,t e ff,L-i, the lemma follows. 

Appendix D 

Derivation of the upper bounds on cr^ K within the 

DIFFERENT INTERVALS 

A. Moderately weak interference regime 

We upper bound er^ K for all values of j3 within each of the 
four intervals. Recall that in the moderately weak interference 
regime SNR~ 1/3 < g 2 < 1. Define 6 = (2c + 8)/ log(SNR), 
where c > is some constant. 

Interval 1 : < |/3| < 1/2 

In this interval the choice a\ — sign(/3) is optimal due 
to ((86>. Therefore, for all \f3\ < 1/2 we have 0~ ai ) 2 > 1/4. 
Thus, 

2 SNR 

^HK > — i— • 
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Interval 2 : 1/2 < |/3| < vIslSNR 



1/4-5/2 



/2 



Define 



Since \/3\ > 1/2 we can express it as /3 = q + <p with 
q £ Z \ and <p £ [-1/2, 1/2). We can further lower bound 

ct hk as 

cthk > min ( ((<P + Q - fl i) 2 + (Q9 - °3 + <pg) 2 ) SNR 
= min ( (if 2 + (qg - a 3 + >pg) 2 ) SNr) . (118) 

The minimizing value of <p (ignoring the constraint ip £ 
[—1/2, 1/2)) is found by derivation to be 



%(q) 



l-2 b 2-2 b - 2fc ' 



q q q 



mod [0,2), (125) 



where I = [—1,1) and the sum in ( 1125b is a Minkowski sum. 
It is easy to verify that 



1 + .9 2 

Substituting ip* into ( 111 8b gives 

2 . / 1 

cr HK > mm 



(qg-a 3 ). 



i,a3 \ 1 + g 2 
1 



{qg-a 3 ) 2 SNR 

> - mini (qg-a 3 ) 2 SNR 

For b = 1,2, ... , [l/31og(SNR)l we define the sets 

g b = {g:g£ [2~ b , 2" 6 + 1 ) } , 
and the quantities 

to ^ V / 2~^SNR 1 / 4 ~ 5 / 2 , 
d> b 4 -J =S NR- 1 / 4 -^ 2 . 



(119) 



(120) 



Let Sb be the set of all values of g £ Qb such that the inequality 

\qg-a 3 \<<5> b (121) 

has at least one solution with < \q\ < q max ,b and a 3 £ Z. 
Note that since q — f3 + ip and we assume in this interval that 

1/2 < |/3| < v /|^|SNR 1/4 - ,5/2 /2, we have 

\q\ < v4^jSNR 1/4 ~ 5/2 . 
Thus, for all g £ Gb and j3 in the considered interval, we have 

\q\ < <?max,b- 



Let Sb = Gb\Sb- Using ( II 19b . we have that for all g £ Sb and 
(3 in the considered interval 

^hk > l^SHR 

1 SNR 172 " 5 

- 2 2- fc +! 

lSNR l/2-«5 

> = — . (122) 



The condition (1121b . which defines the set Sb, can be written 
equivalently as 



\q-2 b g-2 b a 3 \ < 2 b $ b . 



(123) 



Define g = 2 b ■ g, and note that for all j e ft we have 
g £ [1,2). With this notation, (1123b becomes 

2 b a 3 



fJ 



< 2 



b®b 



(124) 



L9m«J 

S b ^2- b |J %(q). 
9=1 



(126) 



Since [^-\ = for all < q < 2 b ~\ for all values of q in 
this range we have 



TM 



mod [0, 2) 



C [2 b $ b l] mod [0, 2) 

= r b (i). 



(127) 



Therefore, 



/ (2 b - 1 -i \ I L*»xJ 

Sb c 2 b |J T b (g) U |J T 6 (g) 

\ V <?=! / \g=2"-l 

= 2- b [76(1) U ( Q 7i(«) | | (128) 

The Lebesgue measure of Sb is bounded by 

n{S b )=Vd(S b ) 

<2- b Vol(T b (l)) + ^ Vol(7i( g )) 

/ LwJ 

< 2~ b 2 • 2 b $ b + ^2 







2 b 





<2$6 + 2$ b 2 2 f^ 



6^ 



< 2$ b + 8 • 2-°$ b g maX)b 
= 2$ b + 8-2" b SNR"' 5 

= V2 ■ 2 b/2 SNR~ 1/4 ~ 5/2 + 8 • 2- b SNR~' 5 . (129) 
We can now upper bound the measure of the outage set 

ri/31og(SNR)l 

s = |J s b , 

6=1 

of all values of SNR~ 1/3 < g < 1 for which ( fT22l does not 
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necessarily hold, as 

ri/31og(SNR)] 



6=1 



l/31og(SNR) + l 
6=1 

l/31og(SNR) + l 



6=1 



Using the identity 



6=1 



P-1 



which is valid for all p ^ 1, and the fact that X^6=i 2 b < 1, 
we have 

fi(S) < V2SNR- 1/4 - s/2 V2SNR 1/6 + 8SNR _<5 

v 2 — 1 



< 7SNR" 5/2 +8SNR-* 

< 16SNR~ 5/2 . 



(130) 



Substituting S = (2c + 8)/ log(SNR) into ( fT22b and ( fT30l ) we 
see that in the interval 1/2 < |/3| < VMSNR 1/4_<5/2 /2 for 
all values of SNR -1 ^ 3 < \g\ < 1 except for an outage set with 
measure not greater than 2~ c we have 

2~ 2c SNR 1/2 



^HK > 



4 -2 5 



Interval 3 : VMSNR 



? l/4-<5/2 



/2 < |/3| < SNR^/VSlfll 



Since SNR 1/3 < g 2 < 1 and we assumed SNR > 4, we 



have 



5 2 SNR- 1 > 



5 2 SNR 



(131) 



Note that CCD) continues to hold for all g 2 > SNR 1/2 . This 
will be useful in the weak interference regime. For all values 
of 1/31 in this interval 



P 



vVSNR-l 



< 



< 



SNR 1 / 4 / y/8\g~\ 
vVSNR-1 
SNR 1 / 4 



<\\9\ 
1 



3/2 SNR -l/4 



and hence, using d85l l, the optimal value of a 2 is 

a 2 = /3/vVSNR- lj =0. 

Therefore, using the fact that S = (2c + 8)/ log(SNR), we can 
upper bound (T84t as 

2 -2c SNR l/2 



^•2 ^ 



/3 2 SNR 
ff 2 SNR- 1 - 4 ■ 2 s 



> 



Interval 4 : SNR 1/4 / A 



\9\ < 



In this interval 



2 > K ~ R 2 > 1SNR 

CT H K > > 7J- 



1/2 



B. Weak Interference Regime 

We upper bound cr 2 ^ for all values of j3 within each of the 



four intervals. Recall that in this regime SNR 
SNR" 1/3 . 

Interval 1 : < \f}\ < 1/2 



-1/2 



< 9* < 



As a\ ^ 0, in this interval (/3 — a\) 

SNR 



> 1/4. Thus, 



'HIC 



> 



Interval 2 : 1/2 < \ f}\ < l/(2\g\) 



In this interval 03 = [fig] 

<J 2 hk > (/3ff) 2 SNR > 



0. Thus, 

5 2 SNR 



Interval 3 : 1/(2|<?|) < \/3\ < yVSNR/8 
Under our regular assumption that SNR > 4, for all values 
of 1/31 in this interval we have 



^5 2 SNR-1 



< 



< 



Vff 2 SNR/8 
V^SNR - T 
y/g 2 SNR/8 
V.9 2 SNR/2 



< 



2' 

where the second inequality follows from (1 1 3 1 i - Thus, the 
optimal choice for a 2 is 



a 2 



/3/V.9 2 SNR- 1 



0. 



(132) 



Therefore, (184b can be lower bounded by 



Interval 4 : vVSNR/8 < 
In this interval 



2 > ^ 2> .9 2 SNR 



Appendix E 
Proof of Theorem[T"21 

For the proof we will need a key result from the field of 
Diophantine approximation which is due to Kleinbock and 
Margulis. The following theorem is a special case of ||57] 
Theorem A]. 

Theorem 14: Let fx, f 2 , ■ ■ ■ , /k be real analytic functions 
in h € U, U a domain in K m , which together with 1 are 
linearly independent over K, and define the manifold 

M = {[h(h) ■■■ f K (h)] : hew}. 
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For almost every hGM and any 5 > the inequality 

max \qhp — aA < \q\~ 7 <~ S 
t=l,...,K 



(133) 



has at most finitely many solutions (q, a) G Z x Z 



A" 



For the proof of Theorem Q~2] we will need a corollary of 
Theorem [14] 

Corollary 7: Let ft, ft, ■ ■ ■ , fx be functions from R m to 
R satisfying the following conditions: 

1) ft for i = 1, . . . , K is analytic in R"\ 

2) 1, /i, . . . , fx are linearly independent over R. 

Let V = {h G R m : /i(h) = 0} and V(e) = V + B(0, e), 
where the sum is a Minkowski sum and B(0, e) is an Tri- 
dimensional closed ball with some radius e > 0. Define the set 
U{e) = R m \X>(e), the set of functions ft(h) = / fc (h)//i(h) 
from U(e) to R for k = 2, . . . , K, and the manifold 

7W(e) = {[/ 2 (h) ••• : he«(e)}. (134) 

For all e > 0, almost every h G M.'{e), and any 8 > the 
inequality 



max 

=1,...,K-1 1 



at 



< 



(135) 



has at most finitely many solutions (q,a) eZx % K . 

Proof of Corollary^ We would like to apply Theorem[T4l 
for the set of functions ft, . . . , fx from U(e) to R. To that end 
we have to show that for all e > the functions ft , . . . , fx are 
analytic in U(e) and together with 1 are linearly independent 
over R. 

The reciprocal of an analytic function that is nowhere zero 
is analytic. Thus, for any e > the function l//i(h) is 
analytic in U(e), In addition the product of two analytic 
functions is analytic. Therefore, for any e > the functions 
ft = A(h) • (l//i(h)) are analytic in U{e) for k = 2, . . . , K. 

We show that the functions 1, ft, ■ ■ ■ , fx from U(e) to R are 
linearly independent for all e > by contradiction. Assume 
they are linearly dependent. Thus, there a exists a measurable 
set S EU(e) and a set of coefficients {ti(e), ...,tjf(e)}£l 
not all zero such that Vh G S 



t^.l+t^^M + .-.+txie)-^ 
/i(h) /i(h) 



0. 



This implies that Vh G S 

• 1 + i x (e) • A(h) + i 2 (e) • / 2 (h) + • • • +t K (e) ■ f K (h) = 0, 

in contradiction to the assumption that the functions 
1, ft, . . . , fx from R m to R are linearly independent over R. 

We can therefore apply Theorem [14] to the set of functions 
fi,---,fK from U(e) to R for all e > 0, and the corollary 
follows. ■ 

We are now ready to prove Theorem Q~2] Define the sets T>, 
T>(e) and U(e) as in Corollary [7] and the manifold 

M(e) = {[ft(h) ■■■ f K (h)] : heU(e)}. 



We begin by showing that for any e > the number of DoF 
offered by the first computation rate is upper bounded by 

1 



-^cornp,! 



< 



K 



(136) 



for almost every h G M.(e). Then we take e to zero in order 
to show that this holds for almost every 

hGTW = {[/i(h) ••• f K (h)] : hel™}. 

Consider the manifold Ai(e) for some e > 0. Note that 
hi — /i(h) 7^ for any h G Ai(e), and we can therefore 
define h = h/hi. We have hi = 1 and [J12 ■■■ hx] G M(e), 
where M. (e) is the manifold from ( 1134b in Corollary [7] 

The channel (O is equivalent to the channel 



1 \ - , 1 

y = 7— y = xi + ) h k x k + — z. 

hi ^— ' /ii 

1 fe^i 1 



(137) 



Let a be a vector of integer coefficients, and j3 be the scaling 
factor used by the receiver in order to compute the lattice 
equation v = EfcLi a fetfe] mod A, see Section [HI] The 
effective noise for computing the equation v with coefficient 
vector a is 



z eff (h,a, (3) = {(i- oi)xi 



y^([3hk - a fe )x fc + — z, 
lii 



and its effective variance is given by 

a e 2 ff (h,a,/3) = (/3- ai ) 2 SNR 



(138) 



Recall that 



-Rcomp,i = max - log 



SNR 



X2"° Va 2 ff (h,a,/3) 
= i log (SNR) - 1 log (^rnin CT 2 ff (h, a, . (139) 

Thus, in order to obtain an upper bound on i? C omp,i we need 
to lower bound cr 2 ff (h, a, 0) for all values of j3 G R and a G 
Z K \0. Let 

h* = max /ifc, 
fe=i,...,if 



and 



k* = arg max h k . 

fe=l,...,K 



Note that if |/3| < l/(2/i*) the minimizing corresponding 
choice of integers a\,...,ax in d 138b is afc« = sign(/3), 
and a k = for all k 7^ k*. This in turn, implies that for 

\P\ < l/(2h*) we have 

SNR 



a 2 ff (h,a,/3) > (rJh* - sign(/?))^SNR > 



(140) 



which means d C omp,i = 0. Thus, in order to obtain a positive 
number of DoF, |/3| must be greater than l/(2h*). 

If l/(2h*) < \f3\ < 1/2, then the minimizing value of 
ai in Cl 1 3 81 > is ai — 0. This implies that for all values of 

l/(2/i*) < \(3\ < 1/2 we have 

SNR 



a e 2 ff (h,a,/3) >/3 2 SNR> 



A(h*y 



(141) 
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which also means d C omp,i = 0. Thus, in order to obtain a 
positive number of DoF, \j3\ must be greater than 1/2. 

Hence, in order to lower bound (11381 in the limit of very 
high SNR, it suffices to limit the optimization space of (3 to 
\(3\ > 1/2. For such values, (3 can be written in the form 
(3 = q + (p where ip e [-1/2, 1/2), and q e Z \ is the 
nearest integer to (3. 

For any \<p\ < 1/2, q e Z \ and a e Z K \0 we have 

offf (h, a, g, 93) 
= (<p + g-0i) 2 SNR 

+ max (qhk - a k iph k ) 2 SNR + ^ ^ 

> ip 2 SNR + max (qh k -a k +wh k ) 2 SNR + 

k^i y r 1 /li 2 



(qh k - a k + <ph k ) 2 ) SNR 



max , ^ 



|2^i| 2 



(142) 

We further bound (11421 by substituting the minimizing value 
of (f for each k ^ 1. It follows by simple differentiation that 
for each fc 7^ 1 the minimum is obtained at 

-h k 



<p*(k) 



(qh k - Ofc). 



Substituting (fc) into ( 11421 i yields 



a eff (h,a, 9 ^)>max Vi + -, 



1 - - -2r„n , <T 



(g/i fc -a fe ) 2 SNR 



> 



max I qh k — a k | SNR 



= c (h)-max|( ? / lfc - afc | 2 SNR+-f^, (143) 
fc^l |2/ii| z 

where co(h) > is some constant independent of the SNR. 

Consider the limit of SNR — > 00, and assume \q\ is upper 
bounded by some finite integer qo > 0. Then, for almost every 
he M.(e) there exists a constant ci(h, go) > 0, independent 
of the SNR, for which 



ma,x\qh k - a k \ > a(h,q ) 



(144) 



for all < \q\ < q and a e Z K ^ 1 . Note that h does not 
satisfy (11441 > only if all elements of h are rational. Substi- 
tuting (I144t into (I1431 > gives <r 2 ff (h, a, g, cp) > c 2 (h, g )SNR 
which means that the DoF is zero. Therefore, in order to get 
a positive DoF, q must tend to infinity when the SNR tends 
to infinity. 

Any positive integer \q\ can be expressed as \q\ = SNR 7 for 
some 7 > 0. From Corollary [7] we know that for any e > 0, 
for almost every he M.(e) and any 5 > 0, for q large enough 
we have 



max\qhk - Ofc| > \q\ K - 1 ' 

k^l 



SNR" — 



-76 



(145) 



Thus, for \q\ large enough and almost every h £ Al(e) we 
have 

a 2 s (h, a, q, tp) > max ^c 2 1 (h)SNR 1 ~" ! ^ T ~ 27 ' 5 , c 3 (h)SNR 27 

(146) 



where 03(h) > is a constant independent of the SNR. 
Minimizing j 146b with respect to 7 gives 

K - 1 



' 2(K + SK — 5) ' 

Hence, for all q E Z, ip e [-1/2, 1/2), a e Z K \0 and almost 
every h 6 M(e) 



cr 2 ff (h, a, g, 93) > c 4 (h)SNR«+ i ^-D 



(147) 



where 04(h) > is also a constant independent of the SNR. 
Substituting into d 1 39b gives 

flcomp,! < k^k^ ■ \ lQ g( SNR ) - \ l0g(C4(h)) 

(148) 

for any 5 > 0. Taking <5 — > 0, it follows that the number of 
DoF the best equation offers is upper bounded by 



Um 



comp.l 



< 



1 



snr^oo i log (1 + SNR) ~ K' 



(149) 



for almost every h £ M.(e). Since this holds for all e > 
we can now take e to zero (note that the bound does not 
depend on e). The set T> has measure zero since fx is analytic 
on K m and is not identically zero (otherwise, the set of 
functions 1, /i, • • • , /if is not linearly independent). Note that 
the measure of T>(e) goes to zero as e — > 0, and furthermore 
T> = D e> o'D(e). Therefore, the claim holds for almost every 
heM. 



Appendix F 
Proof of Theorem[T31 

Consider the reference i-user MAC 

L 

Yref = ^9l*l + Z, 



(150) 



e=i 



where z is AWGN with zero mean and unit variance and all 
users are subject to the power constraint ||x^|| 2 < nSNR. Ap- 
plying Corollary |5]to this channel implies that for almost every 
g € Ai the DoF that each optimal computation rate offers is 
1/L. Let i?j^ f (g, a) be the computation rate corresponding 
to the coefficient vector a over the reference MAC (1150) . 
We now show the computation rate of the same coefficient 
vector i? comp (g, a, B) over the effective MAC (|29l is within 
a constant number of bits from i?^ f mp (g, a). 

For the reference channel (1150) the effective noise variance 
for a given a and (3 is 

<7 2 f (g,a,/3) = SNR||/3g-a|| 2 +/3 2 , 



while for the effective L-user MAC d29l the effective variance 
for the same a and (3 is 

L 

a e 2 ff (g,a,/3,B) = SNR]T(/3 5£ - a e ) 2 b 2 fLi + (3 2 . 

Letting b* = maxti....^ b 2 s e and noting that b* > 1 gives 
a 2 f (g,a,/3) < a e 2 ff (g,a,/3,B) < b*af ei (g, a, (3). 
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Since the above inequalities are valid for any value of /3, in 
particular they hold true for the optimal value of /3 and it 
follows that 

Cmp(S> a ) -\log(b*) < i? comp (g,a,B) < i?^ mp (g,a). 

As b* is independent of the SNR it follows that the DoF 
of each computation rate over the reference and effective 
MACs ( 11501 ) and d29l > offer are equal, and, in particular, this 
is the case for the optimal computation rates, thus the theorem 
follows. 
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